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Abstract: Motivated by the couplings of the dilaton in four-dimensional effective actions, 
we investigate the cosmological consequences of a scalar field coupled both to matter and a 
Maxwell-type vector field. The vector field has a background isotropy-violating component. 
New anisotropic scaling solutions which can be responsible for the matter and dark energy 
dominated epochs are identified and explored. For a large parameter region the universe 
expands almost isotropically. Using that the CMB quadrupole is extremely sensitive to 
shear, we constrain the ratio of the matter coupling to the vector coupling to be less than 
10 -5 . Moreover, we identify a large parameter region, corresponding to a strong vector 
coupling regime, yielding exciting and viable cosmologies close to the ACDM limit. 
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1 Introduction 

Scalar fields play an important role in novel high energy physics models and in theoretical 
cosmology. They provide a simple and natural framework for exploring the possibility 
that the current cosmic acceleration is driven by a dynamical dark energy component — 
quintessence - [1]. If not prevented by an unknown symmetry, these fields are expected to 
couple to other matter fields [2]. The phenomenology of the coupling to dark matter has 
been studied extensively in the literature, see for instance [3] and references therein. 

A natural possibility that seems to be overlooked in models for the late-time universe, 
is the coupling of quintessence to vectors. In fact couplings of scalars to Maxwell- type 



- 1 - 



vector fields are quite common in fundamental theories. For instance, there is a scalar- 
Maxwell coupling in the bosonic sector of the supergravity action. The historically first 
theory that included both the scalar-matter and scalar-Maxwell couplings was the original 
Kaluza-Klein model, almost a century ago. Both of these couplings are also present in 
modern higher dimensional theories such as string theories. It should be stressed that 
typically both couplings are present simultaneously; as a concrete example is the low- 
energy effective gravi-dilaton action [4] which was explored in the context of quintessence 
for example in [5]. So far, however, the gauge kinetic coupling to the Maxwell field has 
mostly been ignored in the context of the late-time cosmology. 

The purpose of this work is to explore the cosmological consequences of these couplings 
for the late-time behaviour of our universe. Although it is natural in the context of higher 
dimensional theories that both couplings are present simultaneously, it turns out that the 
cosmology of this doubly-coupled scenario has not been explored in the literature. We 
shall call the model Doubly Coupled Quintessence (DCQ). Both couplings have, however, 
separately been studied previously although in very different contexts. There is a long 
history for studying scalar- matter coupling going back to the early nineties, see [6-8] for 
early works. We shall refer to the scenario proposed by Amendola in 1999 [9] as Standard 
Coupled Quintessence (SCQ). The main feature of SCQ is a new scaling solution responsi- 
ble for the matter dominated epoch. The scalar-Maxwell coupling, on the other hand, has 
not been studied for the late-time cosmology so far, at least not in the context considered 
here. In this paper we are primarily interested in the back-reaction on the geometry of the 
homogenous vector field. In the context of early universe inflation this has been studied 
extensively recently [10-16]. Sourced by the homogenous vector potential, the inflating at- 
tractor exhibits a small stable anisotropic hair in the expansion rate which, together with 
the anisotropic coupling between variables, has interesting phenomenological consequences 
at the perturbative level [17-27]. In this work we shall show that genuinely new behavior 
arise when both couplings are present simultaneously. The model is phenomenologically 
very rich and exhibits new scaling solutions that can be responsible for the matter domi- 
nated and dark energy dominated epochs. Moreover, we identify a large parameter region 
that (quantitatively) yields cosmologies close to the ACDM limit. We put bounds on the 
model parameters observationally. As we shall see, DCQ can be viewed as a generalization 
of SCQ and the dynamical trajectories of the universe in the latter model represent special 
cases of the more general model. 

In this paper we shall assume that spatial homogeneity holds. Note that a spatially ho- 
mogeneous vector field with a non-vanishing background component picks out a preferred 
direction in the universe and thereby violates isotropy. To consistently study the back reac- 
tion of the vector on geometry we shall assume an axisymmetric Bianchi type I metric. This 
spacetime has homogeneous and flat spatial sections and exhibits a space-like symmetry 
axis which can be aligned with the vector field (one Hubble expansion rate in the direction 
parallel to the vector and one in the plane perpendicular to the vector). It is therefore 
the simplest spacetime consistent with the symmetries of the matter sector. Clearly, our 
model violates isotropy, i.e., three-dimensional rotational invariance, one of the pillars of 
the concordance model of cosmology (ACDM). The cosmic microwave background radia- 
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tion (CMB) provides evidence that the universe is remarkably close to isotropy. Still one 
has the interesting possibility for small deviations from the idealized model. In fact, a shear 
at the one percent level today is consistent with the supernova la data [28]. As we shall see, 
however, the CMB provides a far more effective way to constrain the shear in our model, 
and the upper bound will be significantly lower than what is implied by the supernova la 
data. Since our considered spacetime has the flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric as a special case, we can study solutions which are (arbitrarily) close to 
isotropy and, dynamically, arbitrarily close to ACDM. Interestingly, for a large parameter 
region we shall see that solutions with a small shear are dynamically selected. 

A number of anomalies in observational data hinting towards a violation of rotational 
invariance have been reported by different groups. There are three famous CMB anomalies 
reported in the WMAP, namely the surprisingly low quadrupole [29] (first seen in COBE), 
alignment between low multipoles ("axis of evil") [30, 31], and an hemispherical power 
asymmetry [32-34]. Based on WMAP data a uniform bulk flow of galaxy clusters is 
reported ("dark flow") [35, 36]. Other data analyses that indicate a preferred axis include 
analysis of polarization of electromagnetic radiation propagating over cosmological scales 
[37, 38], and a "handedness" in the orientation of galaxies [39]. It is intriguing that several of 
the reported anomalies seem to point out a common axis in the universe [40, 41]. Although 
the statistical significance is an open issue still being debated [42-44], there is no lack of 
proposed models (or mechanisms) that violates isotropy at late times [41, 45-63]. The 
anisotropy is often associated with dark energy and thus imprints in observables usually 
happens around or after the transition to dark energy domination (z < 0.3). 

The theoretically motivated form of our considered model combined with its rich dy- 
namical behavior clearly distinguishes it from other late-time models. The model is char- 
acterized by three dimensionless parameters that controls the strength of the two couplings 
and the shape of the scalar potential. Leaving the parameters free, we explore the cosmolo- 
gies that arise in the different parameter regions. For a subspace of parameter space the 
cosmology is entirely isotropic and dynamically equivalent to SCQ. In general, however, 
there will be one or more epochs in the cosmic history where the vector is tracking the 
background energy density and consequently the universe expands anisotropically. The 
universe may be anistropic both during the matter and dark energy dominated epochs, 
or possibly only during one of those epochs. To determine the dynamically preferrable 
region of parameter space we use the dynamical system approach and we find that there is 
a well defined parameter region where the shear is sufficiently small to be consistent with 
present observations. Interestingly, in this region the background dynamics is close to the 
ACDM limit. Another exotic possibility that occurs in our model (by tuning the model 
parameters) is that the present state of the universe is the global attractor. This is not 
possible in SCQ (without introducing time dependent couplings) since it is incompatible 
with the presence of a matter dominated epoch. In our model this is achieved due to the 
existence of a new scaling solution which is responsible for the matter dominated epoch. 

In principle the homogenous Maxwell field could be identified with a possible uniform 
component of large-scale cosmic magnetic fields. In our model, however, the coupling to 
the scalar field would correspond to an extremely rapid time variation in the fine structure 
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constant which would be clearly incompatible with constraints from observations. In this 
work we do not identify the coupled Maxwell field with the photons of the standard model 
of particle physics, but assume that it represents a possible sub-dominant dark component 
of the universe which couples only indirectly to matter fields, i.e., via the scalar field and 
the metric. The fact that the present energy density of both dark energy and dark matter 
exceed the energy density of baryons by roughly one order of magnitude, suggests that we 
should keep the possibility open for other undiscovered fields with too little energy or too 
weak couplings to be identified with present data. 

This paper is organized in the following way. In section 2 we introduce the cosmological 
model and derive the field equations and the autonomous system. Section 3 is devoted to 
a general analysis of phase space. Based on this we identify a very interesting parameter 
region which we name the strong vector coupling regime and explore the cosmology of in 
section 4. In section 5 we summarize and discuss our results. Readers who only want a 
quick overview of the main results may jump directly to section 5. 

2 Basics of the model 
2.1 Action 

The doubly coupled model is given by the following Einstein frame action for the metric 
g^ u , scalar field 0, vector field and perfect fluids ^f. 



where g is the determinant of the metric, M p is the Planck mass and R is the Ricci scalar. 
The Lagrangians for the scalar field and vector field are, respectively, 



where the field strength of the vector field is defined by = V^A U — V ' V A^. The Jordan 
frame metric g^ v is related to the Einstein frame metric by the conformal transformation 
9[iv = ^{(fyg^u- Note that the perfect fluid Lagrangian, Cm, describes matter fields (^j) 
which are uncoupled in the Jordan frame. In the Einstein frame, on the other hand, these 
are coupled to the scalar field. Note that apart from the vector part (Ca), the action is 
similar to the case of the standard coupled quintessence (SCQ) model [9]. Our doubly 
coupled model can therefore be viewed as a generalization of SCQ. 

Motivated by dimensional reduction of higher dimensional theories such as string the- 
ories [64] we take the potential of the scalar field and the coupling functions to be expo- 
nentials: 




(2.1) 



£ = --^V^V^-V(0) 
C A = -\f\<P)F, u F^, 



(2.2) 



(2.3) 



V(<j>) 



V e m p , 



(2.4) 
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where A, Qa and Qm are constant parameters characterizing the model. Once the com- 
pactification scheme is specified these can be determined; however, in our model we will 
leave them unspecified and explore the cosmologies that arise in the different parameter 
regions. The simplest (and historically first) model of the form (2.1) is the original Kaluza 
Klein theory, starting with a cosmological constant in 4+1 dimensions. The scalar-Maxwell 
coupling is then a result of the compactification, while the matter coupling arise after a 
conformal transformation to the Einstein frame. 1 Note that when referring to SCQ it is 
understood that we mean Coupled Quintessence with the same exponential type matter 
coupling and scalar potential as considered here, i.e., the case studied in [9]. We consider 
a positive potential, i.e., V$ > 0. We set /o and ho to unity. 2 

Variations of the action with respect to g^, <fi and An, lead to the following equations 
of motion: 

M 2 £^ = T#*> + + T^ M ) , (2.7) 

V„F"" = -«^F"V,A (2.9) 

IVLp 

where E$ is the Einstein tensor, T M = g^T^M and □ = g^V^V u - In the limit C A -> 
the equations agree with the standard formulation of coupled quintessence [9] , while in the 
limit Cm — > we recover the equations of the inflationary model considered in [11, 12]. 
The energy momentum tensor is defined 

and for </> and A^ we get: 

= Q(V0) 2 + V(<f>)\ + V^V^, (2.11) 
T$W = -f^F^p^ _ - A f\mFapF^, (2.12) 

where (V^) 2 = g^ 1 'd^(j)d u (j). To keep the matter sector reasonably simple we shall follow 
most works on coupled quintessence (see [3] and references therein) and use a dust model for 
the coupled dark matter field to be introduced in the next section. This phenomenological 
description is a good approximation to any specific model of dark matter that allows a fluid 
description with negligible pressure, heat-flux and anisotropic stress. It is evident from the 
right hand side of the Klein-Gordon equation (2.8) that there is a coupling of the scalar 
field to matter and to the vector field. From the action we note that the matter field is 
uncoupled in the Jordan frame, i.e.: 

y ftiM) = o. (2.13) 



1 In that case the parameters have the values A = v°, Qa = — V6/2 and Qm = v / 6/6. For a discussion of 
the Kaluza Klein model in the Einstein frame, see [65]. 

2 One can absorb /o into the scalar field by the rescaling <j> — } <f> + <f>o. In the limit of vanishing scalar 
field (0 = 0) we demand g M „ = g M , y implying ho = 1. 
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Next, we reformulate the matter field equations (2.8)-(2.9) in terms of the energy momen- 
tum tensors, and write the matter equation (2.13) in the Einstein frame: 3 

V M T£«0 = _ Q J^A Vv<t> _ Q M ^V v <f>, (2.14) 

V M T^ = +Qa^±V v 4>, (2.15) 
VfiT «M) = +Qm Tm_ Vu ^ (216) 

On this form it is clear that the total energy-momentum tensor is covariantly conserved, 
V a^Tv + + ) = 0. It is also manifest that the couplings lead to exchange 

of energy and momentum both between the scalar field and the vector field, and between 
the scalar field and the matter field. The rate of exchange is determined by the coupling 
constants Qa and Qm, respectively. Note that matter "feels" the vector potential only 
indirectly via quintessence and the metric. The equations of motion (2.14)-(2.16) give 
essentially a fluid picture of the doubly coupled dark energy. 

2.2 The cosmological model 

In this section we shall introduce the cosmological model to be considered in this paper, 
and write down the field equations imposing the simplest possible metric compatible with 
the matter source. In addition to the scalar and vector field, we shall consider the case 
of two matter fields, pressure-free matter (m) and radiation (r). Thus we have = 
j,i4m) _|_ywO_ ^ g ra di a tion has a tracefree energy momentum tensor, radiation is decoupled 
from quintessence. Qm as defined in (2.6) therefore represents a universal coupling to all 
matter fields not being of radiation type. In the context of string theory, however, it is 
expected that quintessence couples differently to different types of matter fields [5]. The 
usual treatment in the literature is to assume that baryons are uncoupled while dark matter 
has one single coupling [3]. We shall neglect baryons in this paper. Since they are both 
uncoupled and subdominant, this should not change our results at a qualitative level, only 
slightly quantitatively. Alternatively one can consider Qm as an extremely small universal 
coupling to both dark matter and baryons (which obviously is a less interesting case since 
there is no hope of confirming such a small coupling observationally) . 

While radiation, pressure-free matter and quintessence are modeled as perfect fluids, 
the homogeneous vector field violates three-dimensional rotational invariance. Still, the 
matter sector possesses a rotational symmetry in the space-like plane orthogonal to the 
vector potential An- We shall consider a spatially flat spacetime with the same symmetries 
as in the matter sector. We choose a coordinate chart such that the symmetry axis is 
aligned with with the x-direction: 

ds 2 = _ dt 2 + e 2a(t) L~Mt) dx 2 + e 2a(t) {dy 2 + ^2)1 _ (2 J?) 



3 Equations (2.14) and (2.15) are obtained by covariant differentiation of (2.11) and (2.12) and using 
the equations of motions (2.8) and (2.9). The identity V 7 F a( 3 + V a Fp 1 + V pF la — and the relations 
Cm = h 4 ((j>)£.M and T£ = h 4 (cf>)f>t are useful. 
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Note that the function e a represents the deviation from isotropy, ie. the shear degree of 
freedom, while e a is the isotropic part. This is a locally rotational symmetric (LRS) Bianchi 
type I spacetime, the simplest geometry consistent with the matter sector. Our considered 
geometry and field configuration, ie. a metric possessing a rotational symmetry in the pure 
space-like plane orthogonal to the vector, is one often assumed in the literature (for an 
incomplete list, see [10, 11, 14-19, 22, 66, 67]). In a previous work [12] we discussed the 
stability and naturalness of this setup. More specifically we considered the most general 
spatially flat and homogenous universe model, which has one more shear degrees of freedom, 
and showed that the expansion normalized shear in the space-like plane orthogonal to the 
vector decays exponentially fast. 4 Thus if the assumed LRS symmetry is broken initially, 
the metric quickly evolves towards an asymptotic limit that agrees with our considered 
setup. 

To quantify the degree of anisotropy in the expansion rate we introduce the dimen- 
sionless shear: 

S = -, (2.18) 

a 

where a "dot" represents differentiation with respect to cosmic time t. Note that an initial 
co-moving sphere will evolve into an ellipsoid that can be either oblate (X > 0) or prolate 
(X < 0) depending on the type of matter content. 5 When the shear is sourced by a uniform 
magnetic or electric field the expansion will be oblate. This is also the case for our coupled 
Maxwell field. 

We shall now introduce a notation that resembles the one used in FLRW cosmology. 
We let a denote the geometric mean of the three scale factors (two of them are equal) : 

a = e a , (2.19) 

and H the mean of the three expansion rates: 

H = a. (2.20) 

Note that H = a/ a. We shall refer to a and H as the mean scale factor and mean expansion 
rate, respectively (although the former actually is the geometric mean). Expressed in terms 
a and H we can define several quantities on the same form as in the FLRW cosmology. 
Note that the proper volume scales proportionally to a 3 , and we define acceleration as 
a > 0. We also define the deceleration parameter in the usual way: 

^-1-^, (2-21) 

and the usual implications q < 4^ a > follows. Next we define a "redshift" time 
parameter on the usual form: 

z = -l + —, (2.22) 
a 



4 This result holds also in a dust or radiation dominated universe, although the decay is quickest during 
an inflationary phase. More concretely, it decays as e _ ' 2_<J ' a where q is the deceleration parameter (2.21). 

The spacetime (2.17) is sometimes referred to as the "ellipsoidal universe". This is, however, a bit 
misleading since the spatial sections are flat, and not "ellipsoidal". The "ellipsoidal universe" can easily be 
interpreted as the Bianchi type IX model which truely has 3-dimensional "ellipsoids" as spatial sections. 
We will thus refrain from using this term here. 
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where ao is the scale factor today. Note that, since the expansion is anisotropic, the 
redshift of photons emitted at a constant time slice (such as the last scattering surface) 
will depend on the direction of the incoming photon. In the appendix we show that, to first 
order in the anisotropy, z is the redshift averaged over all directions, see (C.13) and the 
following comments. In the cosmological interesting case of a small shear, z can therefore 
be interpreted as the average redshift. Thus we have generalized the notion of a "redshift" 
time parameter, z, to the Bianchi I spacetime and demonstrated a clear physical content. 
This allows us to refer to a well defined redshift time-parameter. 

Next, let us consider the matter sector. We assume co-moving fluids and decompose 
the energy momentum tensor relative to a congruence of fundamental observers with four 
velocity in the standard way [68]: 

T^ u = {p + p)u^u u + pg^ + s^, ^ ^ 

Here, the energy density p and pressure p characterize the perfect part, while the imperfect 
part is described by the tensor s^ u . While matter, radiation and quintessence are perfect 
fluids, the vector field is imperfect and contributes both to the isotropic pressure p and 
the anisotropic stress s^v- In this paper we will uniquely work in the coordinate basis 
defined by the metric (2.17), and it is understood that the following tensor components are 
always with respect to this basis. In the coordinate frame the mixed components for our 
considered matter sector can be written: 

T% = diag(-p, p - 2tt, p + it, p + vr), (2.24) 

where it represents the anisotropic stress which is related to the energy density of the vector 
field. We split the energy density and pressure in the contributions from the scalar field 
((f)), vector field (A), non-relativistic matter (m) and radiation (r): 

P = P4> + PA + Pm + Pr, 
P = P<P+PA+Pm+Pr- 

Then, let us consider each of the matter fields. We will assume the standard matter is 
pressure- less, i.e., p m = 0, while radiation satisfies the standard equation of state p r = \p r . 
The energy density and pressure of the scalar field takes the standard form 

P^ = \<P 2 + V(<P), p* = ^ 2 -V(<f>), (2.26) 
with the corresponding equation of state 

For the vector part we shall assume an electric-type field in the x direction (a magnetic 
type field aligned in the same direction is dynamical equivalent to our considered case). In 
the gauge Aq = this corresponds to a vector potential A = A^dx^ = A{t)dx, such that 
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F = (1/2)F /Jill dx^ A dx u = Adt A dx. The energy density, pressure and anisotropic stress of 
the vector field is: 

p A = Ii 2 /V 2a+4CT , PA = l - PA , * = \p A . (2.28) 

Also note that, La = PA for our model. 

Next, we write down the field equations. The gravitational equation (2.7) for our 
metric (2.17) yields 

H + 3H 2 = ^(p-p), (2.30) 
a + 3H& = ^, (2.31) 

while the matter field equations (2.14)-(2.16) become 

Ptf, + SHp^l + Uf) = QAPAjj- - QMPmJj-, (2.32) 

PA + 4(JT + a )pa = —QaPa-^j-, (2.33) 

p m + 3Hp m = QMPmJj-, (2.34) 
p r + AH Pr = 0. (2.35) 

From (2.31) it is clear that the shear is sourced by the anisotropic stress of the vector 
potential. Also note from (2.35) that radiation is decoupled from <j> since the energy- 
momentum tensor is traceless, i.e., g^ u T^ u ^ = 0. From (2.29)-(2.30) it follows that the 
effective (or total) equation of state is 

P 2 H 1 S 2 . 

-p = - 1 ~3lPT^- 2 T^- (2 - 36) 

As described in detail in section 4.4, constraints from the CMB yield |S| <C 1; thus ne- 
glecting £ in (2.36) we define the effective equation of state parameter: 

w e ff = -l-^ = i(2<Z-l). (2.37) 

For acceptable cosmologies we have p ~ co e gp (note that corrections come first to second 
order in £), and we will define the eras dominated by radiation, matter and dark energy by 
^efr — 1/3, uj c g ~ and uj c g < —1/3, respectively. For cosmologies with a large shear, u c g 
has no significance as an equation of state and is merely another deceleration parameter. 
It turns out, however, that for a large parameter region, giving viable cosmology, a small 
shear will be dynamically selected. 
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When deriving equations (2.29)-(2.35) we have neglected the possibility of magnetic 
type components in F^. It turns out, however, that if instead represented a mag- 
netic type field, the only change in the set of equations (2.29)-(2.35) would be the trivial 
substitution Qa — > —Qa- A magnetic type field is therefore dynamically equivalent to our 
considered electric type field. 

2.3 Autonomous system 

To study the phase space structure we now introduce dimensionless variables and write 
down the equations of motion as an autonomous system of equations. The density param- 
eters are defined 

n < s smrf < 2 - 38) 

where the i denote the vector field (A), scalar field ((/>), matter field (m) or radiation (r). 
The Hamiltonian constraint equation (2.29) can then be written on the generic form 

1 = S 2 + fl A + + Mm + ^r, (2.39) 

where X is the shear degree of freedom (2.18). It is useful to split the scalar in the kinetic 
energy and potential: 

= -X 2 + U v , (2.40) 
o 



where 



and 



X=— — (2.41) 

HM p v ; 



We will use the constraint (2.39) to eliminate V(4>) from the equations of motion. The 
system can then be written as a autonomous set of first order differential equations in 
terms of the independent variables: X, S, £Ia, and £l r - It will prove useful to express 
the equation of state parameters, defined in (2.27) and (2.37), in terms of these variables: 

x 2 - m v 

-* = (2 - 43) 



= -1 + \x'' + 2E 2 + -S!j + a m + j!ir, (2.44) 
where in the former it is understood that fiy is a function of the independent variables. 
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We now switch to the dimension less time variable a by use of the identity % = H. 
The autonomous equations can then be written: 

^ = (X + A) (^3(£ 2 - 1) + ^X 2 ^j + 2Xn A + 3(2Q A + A)0 A (2.45) 

3 

+ 3A(fl m + fi r ) + — O m X + 20 r X — 3Qm !!m, 



^ = 20 A (S + 1) + S 



da 
dO A 

da 

d£2 m 

da 
d£l r 

da 



3(£ 2 - 1) + ^X 2 + ^fi m + 20 r 



20 A 



3(S 2 - 1) + -X 2 - Q A X + 1 - 2S + 20 A + -O m + 2Q T 



(2.46) 
(2.47) 



O m [-3 + 6S 2 + X 2 + Q A/ X + 40 A + 3Q m + 4fi r )] , (2.48) 
n r [-4 + 6S 2 + X 2 + 4Q A + 30 m + 4fi r )] . (2.49) 



Since we are considering a positive potential (V((f)) > =>■ Oy > ), it follows from (2.39)- 
(2.40) that the dynamical variables are subject to the constraint 

lx 2 + z 2 + n A + n m + n r < 1. (2.50) 

Since > for all i this implies an upper bound on each variable individually. We shall 
often write Qy > instead of referring to (2.50). Finally we also write down the equation 
of motion for the auxiliary variable, fiy, which is also useful: 6 



dQy 



n v (XX + 3 + 3to cS ) = U v (XX + X 2 + 6£ 2 + AQ A + 3Q rn + 40 r ) . (2.51) 



da 

3 Phasespace analysis 

Equipped with the autonomous system we shall now investigate the phase-space structure 
using a dynamical system approach. In section 3.1 we identify the fix-points of the sys- 
tem and discuss some of their main properties. In 3.2 and 3.3 we study the isotropic and 
anisotropic solutions, respectively, in greater detail; investigate their physical properties, 
classify their stabilities and determine the conditions for existence. The purpose is to iden- 
tify, or rule out, new interesting cosmological scenarios. In 3.4, we characterize parameter 
space in terms of the attractors. Finally, in 3.5, we investigate the generality of the stability 
of the inflationary background solutions. 

As we shall see, the phase space of our model is very rich with 7 isotropic and 6 
anisotropic fix-points which are all physically different. Since the properties of these solu- 
tions depend on a three dimensional parameter space (Qm,Qa,X), a mathematically com- 
plete investigation would include extremely complicated algebraic expressions which are not 
very informative. In this paper we are mostly interested in the cosmology, and therefore, 



6 In simulations, for numerical reasons, it is sometimes useful to replace one of the variables in the 
autonomous system by Qv- In particular this is the case if the initial condition for Qv is (very) small. 
Equation (2.51) can also be used as a consistency check that the autonomous system (2.45)-(2.49) is correct. 
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Table 1. The isotropic fix-points. 



whenever it is desirable and possible, we shall use bounds from observations to simplify the 
discussion. We will thus not emphasize mathematical completeness, and will avoid describ- 
ing the cosmology of models having parameters clearly incompatible with observations, for 
instance the case of a strong matter coupling, \Qm\ 3* 1- Using first order approximations 
will enable us to rewrite most expressions on a simple and informative form applicable to 
extract cosmological results. 

3.1 Overview 

The fix-points of the system are found by setting the left-hand side of the autonomous 
system (2.45)-(2.49) equal to zero and solving the algebraic equations. The stability is 
determined, except in degenerate cases, by linearizing the field equations around the fix- 
points, d ^ = M5X 1 , and evaluating the eigenvalues of the matrix A4. 

We find 7 isotropic fix-points (il)-(i7), and 6 anisotropic fix-points (al)-(a6), see tables 
1 and 2 respectively. In these tables the independent variables (X, S, 0^, f2 m , Q r ) are 
specified, together with cj^, u e g and the auxiliary variable Q^. Especially w e ff characterize 
much of the physics of these solutions. A radiation dominated solution must have uj c q ~ 
1/3, a matter dominated solution u c g ~ 0, while w c fj < —1/3 for an accelerated (dark 
energy dominated) solution. In each of the tables the fix-points are sorted roughly in 
terms if decreasing u c g. The isotropic fix-points (il)-(i7) also exist in, and completely 
characterize, the standard coupled quintessence (SCQ). 7 These well known solutions are 
carefully analyzed elsewhere [3, 9], but since we have two more degrees of freedom (E and 
£Ia), the stability will in general change and needs to be reanalyzed. The phase-space of 
our model is richer than SCQ since we also have the anisotropic fix-points (al)-(a6). The 
scaling solutions (a2), (a3), (a4) and (a5) are genuinely new, while (a6) also exists in the 
inflation model studied in [11, 12, 14]. Apart from (al) (which is of little relevance), non 
of them have previously been considered in the context of late-time cosmology. 

7 When we compare our model to CQ it is understood that we mean CQ with the same type of scalar 
potential and matter coupling as in our model, i.e., exponential functions linear in (/>■ The relation between 
our notation and the parameters j3 and fi in [9] are j3 — —^Qm and fj, = -^A, while the relation to the 
notation in [3] is simply A — > —A. 
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Table 2. The anisotropic fix-points. 



In order to have a self explanatory notation we shall often label the fix-points by 
A, (j>, m and/or r to show which fluids they contain. For instance we shall often write 
(a3)^ mr to remind the reader that (a4) contains vector (A), scalar (0), matter (m) and 
radiation (r). Note that X = 0,a = for all isotropic fix-points (il)-(iT) . Five of the 
fix-points turns out to be so important that we have introduced acronyms beyond this 
labeling. This is RDE(i2), A^MDE(a4), 0MDE(i5), A0DE(a6) and 0DE(i7). The two last 
letters in the acronyms denotes "Dominated Epoch"; for instance "Radiation Dominated 
Epoch (RDE(i2))" and "Vectors-Matter Dominated Epoch (A<£MDE(a4))". This type of 
acronym was introduced by Amendola (in [9]) for the fix-point c/>MDE(i5) which represents 
the main feature of SCQ. We find it natural to introduce the same sort of acronyms in this 
generalized model. 

In the following two sections we shall determine the stability and existence of the fix- 
points. In general the stability and existence depend on the constant parameters A, Qa 
and Qm- Our goal is to get insight into the model by considering different regimes for the 
parameters, and identify those that lead to interesting and viable cosmologies. A given 
fix-point exists in the part of parameter space satisfying the following conditions: X and 
X must be real, (Qa, &m, ^r) must be non-negative and the inequality (2.50) must hold. 8 
This gives in principle six different conditions for existence, but in all cases some of them 
are automatically satisfied while others are related. The independent number of conditions 
on (A, Qm, Qa) is therefore reduced to maximum three. Without loss of generality we 

8 These conditions also imply that > 0. 
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consider A > in the following analysis. 9 

Exact solutions The fix-points correspond to exact power law solutions which are 
straight forward to derive. Prom the matter field equations (2.32)-(2.35) we get: 



Pi* a'*, (3.1) 



where 



3 + 3co^-2Q A X^ + Q M X%^, for i = 0, 

4 + 4S + 2Q A X, for i = A, 
3 — QmX, for i = m, 
4, for i = r. 



(3.2) 



The fact that all fluids scales similarly in the fix-points, implies that ^(Qm , Qa, A) = 
PjiQM, Qa, A) where i and j represents two arbitrary fluids of the considered fix-point. 10 
We therefore write = (3. It is then straight forward to integrate up (2.29) and solve for 
the mean scale factor: 

a oc (3.3) 

3.2 Isotropic solutions 

Here we shall discuss the isotropic fix-points given in table 1. As mentioned above these 
fix-points completely characterize SCQ, and we will mention which labels/names they are 
given in [9] and [3]. Note that the stability will in general be different from SCQ since we 
have two more degrees of freedom (S and £Ia)- 

Fixpoint (ildb)^ This fix-point contains only a pure kinetic scalar field. Thus the fluid 
is "stiff" (w e ff = 1), and therefore dynamically irrelevant (since there is no "stiff" epoch in 
the cosmic history). (il+) and (il— ), respectively, are labeled (e) and (d) in [9], and (bl) 
and (b2) in [3]. Like in SCQ the fix-point is unstable, and exist in the entire parameter 
space. 

Fixpoint RDE(i2) This fix-point contains only pure radiation (Q r = 1) and corresponds 
to the radiation dominated epoch both in ACDM and SCQ. We shall refer to it as the 
"Radiation Dominated Epoch (RDE(i2))", and it corresponds to fix-point (cr) in [9] and 
(e) in [3]. As we shall see in section 4.1, cosmologically viable initial conditions must start 
close to this fix-point also in our model. It is (obviously) a saddle and exists in the entire 
parameter space (\,Qm,Qa)- 



9 Negative A is equivalent to the situation 



10 To avoid possible confusion let us elaborate a bit here. For a fix-point that contains radiation, pi = 4 for 
all fluids contained in that fix-point. For fix-points not containing radiation, pi will be parameter dependent 
and in general of course different from 4. The point is that all fluids that exist in a given fix-point scale 
similarly and thus the corresponding coefficients Pi are equal. 
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Fixpoint (i3), mr This fix-point represents a scaling solution containing a pure kinetic 
scalar, matter and radiation. It is labeled (crm) in [9] and (f) in [3]. It exists in the 
parameter region Q M > 1/2. Apparently, with uj c g = 1/3 it is a candidate for the radiation 
dominated epoch. As we shall see in section 4.1, however, initial conditions starting close 
to this fix-point will not realize a sufficiently stable radiation dominated epoch lasting (as 
a minimum) from big-bang nucleosynthesis (BBN) all the way to redshift z ~ 3000. This 
rules out (i3) as a candidate for the radiation dominated epoch. 11 

Let us also comment on the stability. According to the eigenvalues (A. 3), (i3) is stable 
in the intersection of < o, ^ > 4 and Q M > \. Elsewhere it is a saddle. Also in SCQ 
the stability depends on the parameters and (i3) will either be an attractor or a saddle. 

Fixpoint (i4)^ r This fix-point represents a scaling solution with quintessence (w^ = 1/3) 
and radiation. It is labeled (6r) in [9] and (g) in [3]. It exists in the parameter region 
A > 2 and is stable in the intersection of Qa < and AQm > A. 12 Like (i3) it has uj c g = 1/3 
and is also apparently a candidate for the radiation dominated epoch. As shown in section 
4.1, however, initial conditions starting close to this fix-point cannot provide a sufficiently 
stable matter dominated epoch lasting all the way from z ~ 3000 up to transition to dark 
energy domination at low redshift z ~ 1. Therefore also (i4) is ruled out as a candidate for 
the radiation dominated epoch. 13 

Fixpoint 0MDE(i5) This is the well known scaling solution responsible for the matter 
dominated epoch in SCQ. It is quite routinely referred to as the "(^-Matter Dominated 
Epoch (0MDE)", and is labeled (cm) in [9] and (a) in [3]. It exists in the parameter region 
Q\i < 3/2. However, \Qm\ < 0.13 (which corresponds to < 0.011) is required to satisfy 
CMB constraints [70] . The phenomenology of this fix-point is well understood at the (lin- 
ear) perturbative level. For instance matter fluctuations grows less than in the uncoupled 
case, CMB acoustic peaks are shifted to higher multipoles compared to ACDM (due to 
change in the sound horizon) and low multipoles are tilted (due to the integrated Sachs- 
Wolfe effect) [9]. Regarding the stability, it is interesting to note that that the fourth 
eigenvalue in (A. 5) is positive in the parameter region where the fix-point A0MDE(a4) 
exists. This reflects the fact that in the parameter regions where both 0MDE(i5) and 

"Here we should mention that, although for a different reason, this fix-point is also ruled out as a 
candidate for the radiation dominated epoch in the framework of SCQ [3]. In that case the problem is 
the BBN bounds on primordial quintessence (O^ 1 " 1 < 0.045 [69]) which implies Qf/ > 3.7. Such a large 
matter coupling is incompatible with the existence of a subsequent matter dominated epoch (which must 
be provided by 0MDE(i5) in SCQ). Note that this argument does not hold in the case of our more general 
model since it is possible to realize a matter dominated epoch with such a large matter coupling via the 
new scaling solution A<^MDE(a4). 

12 Note that the real part of the two last eigenvalues in (A. 4) are always negative in the parameter region 
where the fix-point exists. 

13 Although for a different reason, also in the framework of SCQ this fix-point is ruled out as candidate 
for the radiation dominated epoch [3]. In that case the problem is that the BBN bounds on a primordial 
quintessence (mentioned in the above paragraph) implies A 2 > 88.9 which is incompatible with a dark 
energy dominated epoch (which must be realized by 0DE(i7)). Note that this argument does not hold in 
the framework of our more general model since a viable dark energy dominated epoch can be realized by 
the scaling solution A</)DE(a6) as long as the the fraction |A/Qa| is small. 
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A<^>MDE(a4) exist, the latter one will be dynamically selected. In section 4.2 we shall see 
that in general A0MDE(a4) will be responsible for the matter dominated epoch, while 
0MDE(i5) will be responsible only in the parameter region where A0MDE(a4) does not 
exist. For \Qm\ < 0.13 t/>MDE(i5) will in general be a saddle followed by a late time ac- 
celerated solution. However, in the parameter region where A0MDE(a4) does not exist, it 

3+2Q 2 

will be a stable attractor provided the additional requirements Qm > and A > 2Qm ' • 
Obviously, this parameter region (notice that A > when Qm <C 1) can be ruled out since 
it will be impossible to realize the late-time accelerated epoch. 

Fixpoint (16). This is a scaling solution with matter and quintessence. It is labeled 
(6m) in [9] and (d) in [3]. In principle (i6) is a candidate for both the matter and dark 
energy dominated epoch. In the framework of SCQ, (i6) can be ruled as a viable candidate 
for any of these periods [3, 9]. Below we show that in our more general model the same 
conclusion holds only for the matter dominated epoch. Quite interestingly, it is possible 
to realize a viable cosmology where (i6) represents an accelerated late-time attractor. The 
fix-point exists in the parameter region where 2Qm{\ — Qm) < 3 and A(A — Qm) > 3, where 
the first condition comes from VLy > and the latter from {l m > 0. 

Then, let us show why (i6) is not a viable candidate for the matter dominated epoch. 
In the limit A 3> \Qm\ an d A > 1 the fix-point is apparently a candidate for the matter 
dominated epoch with uj c r ~ and Q m ~ 1. In order to have a subsequent accelerated 
epoch we must further require 6Qa ^ A so that we have saddle stability (follows from 
a detailed study of the eigenvalues (A. 6) combined with the existence conditions). Since 
A > 1, the only possibility for a subsequent accelerated epoch is provided by A0DE(a6). 
Then, \Qa\ 3> A is required in order to have a small shear in the accelerated era. Combined 
with the above conditions we then have —Qa >A>1. In that case A(/>DE(a6) does not 
exist. Thus it is impossible to realize a subsequent accelerated epoch. This rules out {jft)^ 
as a possible candidate for the matter dominated epoch. 

Next, let us consider (a6) as a candidate for the accelerated epoch. Acceleration 
(w c ff £ [ — 1) — \)) is realized in the parameter region Qm < — ^A. 14 Under this condition 
the two first eigenvalues in (A. 6) is negative. If we additionally have Qa < X ~^9 M and 
Qm < A — j, also the other eigenvalues are negative. In that parameter region (iG)^ 
represents an accelerated scaling solution with matter and quintessence. This includes the 
special case (Qm, A) = (—3.3, 2.2) in which case ~ 0.7 and U m ~ 0.3. This opens up 
the exotic possibility that the present state of the universe represents a global attractor. 
Interestingly, in this scenario the universe has already reached the final attractor solution 
and is not in the middle of the transition between two cosmic epochs. As mentioned above, 
in the framework of SCQ this possibility is inconsistent with the presence of a matter 
dominated epoch since </>MDE(i5) does not exist in the considered point in parameter space 
[3, 9]. In our more general model, however, the matter dominated epoch can instead be 
realized by the scaling solution A0MDE(a4), and, quite remarkably, it is indeed possible 
to realize a viable cosmology on the background level. We consider the details of this 
possibility in section 4.3. 

14 For Qm > A we have u e g < — 1 and the fix-point does not exist. 
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Fixpoint 0DE(i7) This is the well known quintessence attractor which exists when A 2 < 
6 and accelerates when A 2 < 2. It is labeled (a) in [9] and (c) in [3]. Note that it represents 
the late-time accelerated solution both in SCQ and in "conventional" quintessence without 
matter coupling (with exponential type potential). In our more general model it is one 
of three possible late-time attractors that can realize acceleration (together with (16)^ 
and A0DE(a6)). From the eigenvalues (A. 7) it follows that the fix-point is stable and an 
attractor in the intersection of A 2 < 4, Qm > A — f and Qa < ^rr - - 



3.3 Anisotropic solutions 

Here we shall discuss the anisotropic fix-points given in table 2. As mentioned above, (al) 
is a well-known solution of little relevance; (a2),(a3),(a4) and (a5) are genuinely new scaling 
solutions; while (a6) is explored in the context of inflation [11, 12, 14]. Apart from (al), 
non of them are considered in the context of late-time cosmology. 

Fixpoint (al±),* This is a decelerating fix-point containing only a pure kinetic scalar 
field (thus the fluid is "stiff" p = p). Since there is no "stiff" epoch in the cosmic history, 
this solution is dynamically irrelevant. The solution is part of a broader solution commonly 
referred to as Jacobs disc [71]. See also [72] for a discussion of such solutions. Note that 
(e) is a curve of fix-points, satisfying ^X 2 + T? = 1. Notice from the eigenvalues (A. 8) 
that the fix-point is unstable. 

Fixpoint (a2)^ r This is a genuinely new scaling solution containing radiation, vector 
and quintessence. With w c g = 1/3 it is apparently a candidate for the radiation dominated 
epoch. However, in (4.1) we shall see that it is impossible to realize a sufficiently stable 
radiation dominated epoch lasting (as a minimum) from big-bang nucleosynthesis (7bb n ~ 
lMeV) all the way to redshift z ~ 3000. 

The fix-point exists in the parameter region where Qa > and A > \Qa~\-\\J + 16- 
The former follows from Qa > 0, while the latter from $l r > 0. Notice that the fix-point 
does not exist in the parameter region A < 2. The fix-point is a saddle in the major frac- 
tion of the parameter space (Qm, Qa, A), but note that there are certain spots where it is 
a stable attractor, see eigenvalues (A. 9). 

Fixpoint (a3)^^ mr This is a scaling solution with radiation, matter, vector and the 
kinetic part of the scalar. With w e g = 1/3 it is apperantly a candidate for the radiation 
dominated era. However, in section 4.1 we show that it is impossible to realize a subsequent 
viable matter dominated epoch. Thus we can rule out (a3) as a candidate for the radiation 
era. 

The fix-point exists in the intersection of > and ^Q M — QmQa — ?>Q\ > 2. The 
former follows from Qa > 0, while the latter from f2 r > 0. Thus the fix-point does not 
exist for a small matter coupling \Qm\ 1, more precisely it does not exist in the region 
Qm £ (— "^j"^)- Like (a2) this fix-point is in general a saddle, but in certain spots of 
parameter space (Qm,Qa,^) it is a stable attractor, see eigenvalues (A. 10). 
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Fixpoint A0MDE(a4) This is a genuinely new scaling solution with matter, vector and 
the kinetic part of the scalar field. As we shall show below this fix-point is a candidate for 
the matter dominated epoch in the limit of a strong vector coupling. In chapter 4.2 we shall 
see that, whenever it exists, A0MDE(a4) will be the dynamically selected solution during 
the matter dominated epoch. We shall refer to (a4) as the "Vector-</>-Matter Dominated 
Epoch (Ac/>MDE(a4))" . The exact eigenvalues are given in (A. 11). In the analysis below 
we shall focus on the parameter region where A0MDE(a4) represents a matter dominated 
solution. 

The conditions for existence follows from Qa > and O m > 0. For \Qm\ < this 
gives the single condition AQmQa + 2Q m > 1- For a small matter coupling (\Qm\ "C 1) this 
condition reads \Qa\ S> 1 (with the sign of Qa equal to the sign of Qm)- In this region we 
have to lowest order Q m = 1 and S7 c ff = 0, which means that the fix-point is a candidate 
for the matter-dominated epoch. Under the assumption \Qa\ 3> \Qm\ and \Qa\ » 1 we 
expand the characteristics of the fix-point to first order in the small quantities ^p- and ^-: 



X - 1 v — 2( ^ M _ 1 q _ Qm 
2Q a ' ~ 3Q a 6Qi' A ~2Q A ~ 

n i ® M 1 n n o 1 1 Q M 

nm - l ~WA + ^ I A ~ 24Qf ' U + = 1 ' ^~Wa 



?a 



(3.4) 



Thus, to zeroth order, Ac/>MDE(a4) looks like a standard matter dominated epoch with 
£l m = 1. In section 4.4 we shall constraint the ratio Qm/Qa observationally by requiring a 
small S compatible with observations. Applying equations (3.1)-(3.3) we find, in the same 
approximation, that the energy density and the mean scale factor evolves as: 

poca " {L+ sQa\ a<xt 3( 6 Qa>. (3.5) 

Note that the A0MDE(a4) is dynamically equivalent to the matter dominated epoch of 
ACDM in the limit Qm/Qa — > in which case £l m = 1, p oc a~ 3 and a oc t 2 ^ 3 . In the same 
approximation as above the eigenvalues are: 

3 Qm Qm 3 Qm 3 Qm » Qm A \ 

"2 + 4Q A ' IQa 4 + 8Qa' 4 + 8Q a ' 2Qa 2Qa)' l ' J 

The first four eigenvalues are negative under the assumptions of the approximation. The 
last one is positive unless X/Qa ^ 6, i.e., A > 1. Thus we must require X/Qa < 6 to secure 
that A0MDE(a4) is a saddle followed by an accelerated epoch. 

Fixpoint (a5)^ m This is a genuinely new type of scaling solution with matter, vector 
and quintessence. It exists in the intersection of £Ia > 0, Vt m > and Sly ^ 0, which 
yields complicated algebraic expressions in terms of (Qm,Qa, X). In principle it could 
provide interesting cosmological scenarios. Among the possibilities are accelerating scaling 
solutions, an "almost de Sitter" solution or alternatives to the matter dominated epoch. 
In appendix B we have explored the cosmological viability of all these possibilities. It is 
clear from the expression for E in table 2 that such alternatives would require fine-tuning 
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in the parameters to secure a small shear (|S| <C 1) compatible with observations. In 
the appendix we show that even with an appropriate fine tuning, these alternatives are 
not cosmologically viable and can be ruled out. We conclude that (a5) can neither be 
responsible for the matter nor the dark energy dominated epoch. 

Fixpoint Ac/>DE(a6) This is a scaling solution with vector and quintessence that can 
realize the late-time acceleration. It was first discovered in the context of inflation recently 
[11, 14]. Providing a counter-example to the cosmic no-hair theorem it has already gen- 
erated much enthusiasm and attention. In the context of our model it plays the role as a 
possible future attractor together with 0DE(i7). In section 4.2 we shall see that in parame- 
ter regions where both 0DE(i7) and A^MDE(a4) exists, the latter one will be dynamically 
selected. 

The existence of this fix-point relies on the following two conditions: Q,a > and 
Qy > 0. These two conditions can be written A 2 + 2AQ,4 — 4 > and | < x*+2\qI-2 ' 
To avoid a large shear incompatible with observation we must require |E| <C 1 which is 
satisfied in the parameter region Qa 3> A. For this case the existence condition simplifies 
to 2 < XQa- Thus there is a large parameter region where (a6) exist and can be responsible 
for a small shear. 

Under the assumption Q\ 3> 1 and Qa > A we expand the characteristics of the 
fix-point to first order in the small quantities and -q^: 



Qa' 3Qa 3Q|' A 2Q A Qf 

1 A 1 2A 1 

\l m = 0, U r =U, Ufj, = 1 - -— + — J, U e ft = -1 + — — + 8— 2- 

2Q A Qi 3Qa Q% 



(3.7) 



Thus, to zeroth order, A0DE(a6) looks like a de Sitter solution with = 1 and cj e ff = — 1- 
In section 4.4 we shall constraint the ratio X/Qa observationally by requiring a small 
shear compatible with observations. Applying equations (3.1)-(3.3) we find, in the same 
approximation, that the energy density and the mean scale factor evolves as: 

2A Q A 

pot a q a , a oc i * . (3.8) 
In the same approximation the eigenvalues read: 

3 + — 4 + 2— - + IA ^ + IA 3 + 2— 2—^) (3 9) 

Qa ' Qa ' 2 2 Qa ' 2 2 Qa ' Qa Qa J ' 

which shows that under the assumptions of the approximation, and the weak additional 
requirement that \Qm\ is not extremely large, A0DE(a6) is a stable attractor. In other 
words, A0DE(a6) is the global attractor in the parameter region where it realizes a small 
shear compatible with observations. 

3.4 Parameter space 

In this short section we shall use the results derived above to characterize parameter space 
in terms of the attractors. Figure 1 shows the parameter space for the slices Qm = 0.1 
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a)Q_M=0.1 b)Q_M=-3.3 




Figure 1. Figure a) shows the parameter space (\,Qa) of the model for the slice Qm = 0.1. Each 
region is labeled by the fix-point that is an attractor in that region. The green shaded region is 
the accelerated part of the attractors. For Qm = —0.1 the plot is similar except that (i6) is the 
stable fix-point also in the regions labeled (i5) and (a4). Figure b) shows the the same for the slice 
Qm = —3.3, a case to be discussed in section 4.3. 

and Qm = —3.3. Parameter space is divided into different regions each labeled by the 
fix-point that is an attractor in that region. Neglecting the boundaries we note that each 
point in parameter space has one and only one stable attractor. The green shaded region 
represents the part of parameter space where the attractor is accelerated, i.e., w c g < —1/3. 
Obviously, since we have ruled out the possibility of an accelerated saddle, the white region 
is ruled out since it cannot realize a dark energy dominated era. In the next section we 
shall derive further, and stronger, bounds on the parameters by requiring consistency with 
the (observed) CMB quadrupole. For Qm = —0.1 the plot is similar as the case Qm = 0.1 
in figure 1 except that (i6) is the stable fix-point also in the regions labeled (i5) and (a4). 
Note that for the slice Qm = —3.3 the matter-quintessence scaling solution (i6) can realize 
an accelerated attractor. We will discuss this case further in section 4.3. 

Note that we have not considered the boundaries between the different regions of 
parameter space. For these special points a more careful analysis is needed (see, e.g., [73]). 
Furthermore, there may also be regions which contain other, more exotic, asymptotes, 
e.g., attracting closed curves. Although such possibilities in parameter space may pose 
mathematically interesting dynamical solutions we shall not attempt to identify them here 
as they merely represent special cases of our free parameters. 15 

3.5 Generality of the stability of the inflationary attractors 

In spite of the fairly complicated nature of the general system of equations, it is still 
possible to consider - at least heuristically - the stability of the exact solutions with 
respect to generic inhomogeneous perturbations. In order to achieve this we will adopt 

For more complicated Bianchi models such "exotic" behaviour is expected since such behaviours have 
been seen even in ordinary general relativity with perfect fluid [74, 75]. 
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the formalism for studying the dynamics of inhomogeneous cosmologies (Go models) by 
employing expansion-normalised scale-invariant variables developed in [76] . This formalism 
uses orthonormal frames to write the Einstein field equations as an autonomous system 
of evolution equations and constraints. This formalism was successfully used to study the 
initial singuarity in [76], and the instability of inhomogeneous perturbations in the more 
complicated Bianchi type VII^ plane waves [77]. 

We should point out that although the equations in [76] are for a 7-law perfect fluid, 
some of the equations are purely geometrical in nature, i.e., stemming from equations like 
the Jacobi identity. The heuristic argument used here only depends on the background 
solutions and equations of geometric origin, thus the actual form of the matter source will 
be supressed. When that is said, the full set equations for the scenario we are considering 
here is fairly straight-forward to derive by generalising the equations in [76] (see e.g., [78]). 

The orthonormal frame {eo,e a } can be expressed in local coordinates as 

e = ^(d T -N i d i ), e a = ejdi. (3.10) 

The Hubble-normalised frame is defined in general as follows: 

do = -^e , d a = ^e a . (3.11) 

In addition to the deceleration parameter q, it is also necssary to introduce the spatial 
Hubble gradient r a : 

q+l = -^d H, r a = -^d a H. (3.12) 
Introducing the separable volume gauge: 

N = H~ 1 , N* = 0, U a = r a , (3.13) 

enables us to write 

d = d T , d a = Ejd h (3.14) 

where EJ = e^/H. 

The Hubble-normalised state vector for Go cosmologies can then be given by 

X= [Ej, r Q , E a/3 , A a , N af} , Q, P, U a p, Q a f , (3.15) 

where S Q( g is the expansion-normalised shear, A a and N a p are the connection variables, and 
$7, P, H a p, and Q a are matter variables (the total expansion-normalised energy density, 
isotropic pressure, anisotropic stress and energy flux, respectively). The evolution equations 
and the constraints can now be written down in terms of X (see [76]). The equations can 
be written on the form 

d r X = F(X, d;X, didjX), (3.16) 
= C(X,$X), (3.17) 
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where eqs.(3.16) are the evolution equations and eqs. (3.17) are the constraints. 

The frame variables, EJ decouple in the spatially homogeneous subset and thus have 
no direct consequence for the dynamics in this subset. However, for general inhomogeneous 
cosmologies they may have an impact on the dynamics. We are interested in perturbing the 
fix-points found earlier, in particular, the stable fix-points (i7) and (a6) are of particular 
interest. Consider these fixed points, for which we can choose coordinates such that the 
frame variables are: 

Ej = diag(£f , El El). 
The evolution equations and the commutator relations are 

d T Ej = (qdi - ?P a + e^E?) Ej (3.18) 
= 2 (d [a - r [a - A [a ) Eft - e a p 5 N s "iEj, (3.19) 

In a Fermi-propagated frame, W = 0, and for the fix-points q and = diag(— 2S, E, S) 
are constant; thus the evolution equations can be integrated to give: 

Ef oc e (9+ 2S ) T , E\ oc e ( ^ s)r , El oc e (9 ~ E)r . (3.20) 

When perturbing the fix-point solutions with respect to the generic inhomogeneous 
perturbations, the linearised equations will depend on the scale-invariant partial derivatives 
Q X = E l a diX, where E l a are the solutions above. In particular, if / = f(r,x,y,z), then 

d x f oc e^ +2 ^ T d x f, d 2 foce^ T d y f, d 3 foce^ T d z f, (3.21) 

to lowest order. 

To investigate the (in)stability of these spacetimes with respect to general inhomo- 
geneous perturbations we linearize the equations of motion around the exact fix-point 
solutions (assuming that terms like (X — Xo), c? Q X, d a dpX, etc. are small). The inter- 
esting modes stem from the Hubble gradients r a , so to understand the dynamics of the 
perturbations we must linearise these equations explicitly. The evolution equation for r a 
and the vorticity constraint read 

d T r a = (qdP a - tf a + e P ai W) rp + d a q, (3.22) 







£ 



a ^-^)-Hr T (3.23) 



The deceleration parameter, q, is given by 

q = 2T? + 1(0 + 3P) + 2 -A a r a - ^d a r a . 

Lastly, we also need to consider the stability of the curvature variables (which can be 
interpreted as spatially homogeneous modes): 

d T A a = (q5 a p - £ a £ + e a lP K 1 )A^ + \dp(Y> ap + e a ^W), (3.24) 
d T N a P = {q5 {a 5 + 2S ( ^ + 2e^ a W)N^ 5 - 7 ( e T*( a E^ - F^R® + 5 af} W). (3.25) 
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Starting with the latter two equations, we note in the absence of inhomogeneous modes, 
that the fix-points considered here are stable to spatially homogeneous modes as long as: 

(g-4S)<0, (g + 2S)<0. (3.26) 

The inhomogeneous modes are a little bit more difficult to treat since there may be solu- 
tions for which <9jX might diverge even if E % a — > 0. We therefore need to make an assump- 
tion about the inhomogeneties. The inhomogeneous modes enter the evolution equations 
through c? a X = E l a diX. Let us compare with inhomogeneous modes _fk = exp(ik • x), 
where k = (k x , k y , k z ) are in general time dependent: 

d a h = i{Kh)f v . 

Whether such modes grow or decay depends on the term E l a ki. If 

(g + 2£)<0, ( 9 -E)<0, 

then the E l a will decay and, if ki does not grow too quickly, then E l a ki — > 0. Moreover, are 
these conditions satisfied, one can see from the equations above that as r — > oo: 

(& a ,r a ,d a X) -> (0,0,0), 

(in line with [76], this is the condition for future asymptotic silence). In particular, ki grow- 
ing can be interpreted as the inhomogeneities "collapsing" (in terms of the coordinates). 
However, they only have cosmological significance if {E % a ki) do not remain small. This 
implies that inhomogeneties which have modes ki not growing too rapidly will experience a 
" homogeneisation" , or smoothing out, during an inflationary phase for which (q + 2T,) < 
and (q — E)< 0. 

For our most interesting solutions, the fix-points (i7) and (a6), we have (q, E) ~ (—1, 0) 
and (q, S) « (—1, 10~ 5 ), and thus all of the bounds are fulfilled. Consequently, these are 
stable against homogeneous pertubations, as well as inhomogeneous perturbations provided 
that the modes ki do not grow too quickly. 

We should also point out that the equations used here are purely geometric (aside from 
the actual expression for q), thus the heuristic perturbation argument is model independent 
and fairly robust. Thus, we expect similar stability for other models with similar behaviour. 

4 Cosmology of the strong vector coupling regime 

Based on the dynamical system analysis above we shall now identify a cosmologically very 
interesting region of parameter space. This is the strong vector coupling regime defined by 

\Qa\ > 1, \Qa\ > IQwI, \Qa\ > A. (4.1) 

The goal of this section is to understand the cosmology of this parameter region and 
put bounds on the parameters of the model. As we shall see, the strong vector coupling 
regime yields viable cosmologies with important differences compared to Standard Coupled 
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Quintessence (SCQ). Given appropriate initial conditions in the radiation dominated era 
(u c g ~ 1/3), a subsequent matter dominated epoch (w e g ~ 0) will be realized before 
entering the accelerated era (w e fT < — 1/3)- In general the dynamically selected solutions 
in the matter and dark energy dominated epochs will have a small expansion anisotropy. 
The associated shear variable is automatically small in the strong vector coupling regime. 
Consistency with the CMB quadrupole will provide upper bounds on the fractions \Qm/Qa\ 
and |A/Qa|- No additional fine-tuning in any of the parameters is needed to avoid a 
dominating and observationally unacceptable shear. 

As we shall see, the strong vector coupling regime is characterized by only a handful of 
fix-points. These are RDE(i2), A0MDE(a4), 0MDE(i5), (ie)^, A<£DE(a6) and </>DE(i7) 
which we analyzed individually in the previous section. In this section we shall analyze 
the cosmological trajectories that connect these solutions. We shall see that the cosmology 
can easily be understood qualitatively in terms of these fix-points. 

In section 4.1 we identify the initial conditions yielding viable cosmologies. In 4.2 we 
study the possible trajectories starting with viable initial conditions. In 4.3 we consider 
a special case where the present universe, represented by (16)^, is the global attractor. 
In 4.4 we consider the CMB quadrupole and put bounds on the parameters of the model. 
Finally, in 4.5, we consider the sound horizon. 

4.1 Initial conditions 

Among the fix-point solutions in our model several have u c g = 1/3 and are therefore 
natural candidates for the radiation dominated era. These are the "concordance" radiation 
saddle RDE(i2), which has Q r = 1, and the scaling solutions (i3)^ mr , (i4)^ r , (a2)^, , and 
(s^ASmr- In this section we shall argue that the scaling solutions can be ruled out, and 
that cosmologically viable initial conditions must be represented by the concordance saddle 
RDE(i2). We show this by assuming that the matter dominated epoch is realized by either 
A0MDE(a4) or </>MDE(i5). 16 

Viable initial conditions must provide a matter dominated era squeezed in between 
the radiation dominated era and the accelerated epoch. The matter dominated epoch 
must last several e-folds such that matter perturbations can grow and eventually collapse 
non-linearly on scales up to approximately 1 Mpc. Obviously, this puts strong bounds on 
the initial conditions. First note that the matter dominated fix-points A^>MDE(a4) and 
0MDE(i5) have = 1, i.e., the kinetic part of (ft is responsible for the scalar field energy. 
From (2.51) it follows that the potential energy of the scalar field will scale approximately 
as fly ° 3 during a matter dominated era (assuming A < 1). Thus, in order to have 
a matter dominated era lasting several e-folds, Sly must be extremely small at the time 
of matter-radiation equivalence. This observation implies that a viable candidate for the 
radiation dominated era must have uj^ = 1. Otherwise, i.e., if oj^ < 1, the Cly "handed 
over" from the radiation dominated era will be way too large (for any reasonable parameters 
\QmiQa) to provide a matter dominated epoch sufficiently stable to last from z ~ 3000 

16 In section 3 we ruled out (i6) , and (a5) A , as candidates for the matter dominated epoch. 
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until z ~ 1. We can therefore rule out (i4)0 r and (a2) A<7ir , as candidates for the radiation 
dominated era. 

Next, let us consider the scaling solution (i3)0 mr which has 00^ = 1 and Q^ = 1/{QQ M ). 
The bound on Q^ at big bang nucleosynthesis is [69]: 

f^ bn < 0.045, (4.2) 

which means that Qlj > 3.7 given that (i3) is responsible for the radiation era. In standard 
coupled quintessence such a large \Qm\ can be ruled out since it is impossible to realize 
the subsequent matter dominated era via 0MDE(i5) which only exists for Q M < 3/2. 
Apparently, our more general model can realize a subsequent matter dominated era via 
A0MDE(a4) with a large \ Qm\- It turns out, however, to be impossible is to get a sufficiently 
stable radiation dominated era that lasts, as a minimum, from big bang nucleosynthesis all 
the way to redshift z ~ 3000 (where radiation-matter equivalence is supposed to occur). 
To see this, first note that \Qa\ ^ \Qm\ is required to have a small shear (£ ~ Qm/Qa) 
during A^>MDE(a4). Furthermore Qa and Qm must have the same sign for A0MDE(a4) 
to exist, i.e., Qa/Qm 3* 1- In section 4.4 we shall see that consistency with the CMB 
quadrupole requires Qa/Qm ^ 10 5 . Next, consider the field equation for Qa, equation 
(2.47), that includes a term — QaX on the right hand side. Evaluated at (i3) this term is 
Qa/Qm- Such a large term would cause the vector to grow absurdly fast, Qa tx a 10 °, at 
least initially when it is subdominant. After the initial rapid grow, simulations shows that 
Qa will start to oscillate and in fact stay subdominant. This will, however, trigger Q m 
to start growing and the radiation era is replaced by the matter era within a few e- folds. 
To have a radiation era stable from BBN to matter-radiation equality would require an 
unreasonable fine tuning in the initial condition for Qa- The only reasonable possibility 
is to put Qa = as initial condition, but in that case A0MDE(a4) cannot be realized 
and consequently there is no matter dominated epoch. Clearly, we can rule out (i3) as a 
candidate for the radiation dominated era. 

Finally, let us consider the scaling solution (aS)^^ ■ From the conditions for existence 
given in section 3.2, it follows that the fix-point does not exist in the strong vector coupling 
regime defined in (4.1). Even outside the strong vector coupling regime (a3) can be ruled 
out since the existence of t/>MDE(i5) is incompatible with the BBN bound (4.2) and the 
shear of A0MDE(a4) is large in the parameter region where both (a3) and (a4) exists. 

Using only weak assumptions, we have showed that initial conditions must be chosen 
close to the concordance radiation saddle RDE(i2) to get a viable cosmology. It follows that 
the potential energy of the scalar field will scale approximately as a cosmological constant 
during the radiation era (Qy oc a 4 ), and in fact also during the matter era (Qy oc a 3 ). It 
follows that Qy must be extremely small initially. Exactly the same sort of fine-tuning is 
required also in other quintessence models with the exponential-type potential. 

4.2 Dynamics 

Although our model is phenomenologically very rich, we shall now see that in the strong 
vector coupling regime the dynamics is characterized by only a handful of fix-points and can 
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Table 3. The fix-points characterizing the strong vector coupling regime. Several quantities are 
approximated using the assumptions of the strong vector coupling regime (4.1). 



more easily be understood qualitatively. In addition to the constraints on the parameters 
coming from the strong vector coupling regime, see (4.1), we shall assume that \Qm\ an d 
A are smaller than unity in this section. Phase space is then characterized by the fix-point 
solutions RDE(i2), A(/>MDE(a4), </>MDE(i5) , A0DE(a6) and 0DE(i7). Properties of these 
fix-points are summarized in table 3. 

We should emphasize that although we consider A and \Qm\ smaller than unity in 
this section, there are also large parameter regions yielding viable cosmologies when they 
are larger than unity. For A and/or \Qm\ larger than unity, the situation is a bit more 
complicated, however, since one must carefully avoid the problematic regions of parameter 
space discussed in section 3, e.g., regions where 0MDE(i5), (a5)^ m or (i6)^ m represent 
matter dominated attractors (none of these undesirable situations exist when | Qm | and A 
are smaller than unity). Another difference is that 0MDE(i5) and 0DE(i7) do not exist 
when Qm and A are large (see table 3) . An interesting possibility, to be discussed in section 
4.3, that occurs when \Qm\ and A are larger than unity, is that (i6L m represents the present 
universe. Apart from special cases involving (i6), the following discussion covers the main 
features of the cosmology of the strong vector coupling regime. Several of the scenarios 
discussed in this section also exist when \Qm\ or A are larger than unity. 

Viable solutions must realize a matter dominated era squeezed in between the radiation 
and dark energy dominated epochs; in other words the following sequence must be realized: 



(^cff-1/3) -> (w cff ~0) -> (u cS < -1/3). 
In our model this sequence can be achieved in various ways: 

/ A<£MDE(a4) \ / A<£DE(a6) 



(RDE(i2) 



or 



or 



(4.3) 



(4.4) 



\ </>DE(i7) 



\ </>MDE(i5) 

As shown schematically, the matter dominated epoch can be realized either by the anisotropic 
scaling solution A</>MDE(a4) or the isotropic scaling solution 0MDE(i5), while the dark en- 
ergy dominated epoch can be realized either by the anisotropic scaling solution A^>DE(a6) 
or the isotropic solution 0DE(i7). This yields four qualitatively different trajectories. Ex- 
amples of each of these trajectories are given in figure 2. In general the anisotropic solutions 
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are dynamically preferred over the isotropic ones. If both A^MDE(a4) and </>MDE(i5) ex- 
ists, the former will be dynamically selected in the matter dominated epoch. Similarly, if 
both A0DE(a6) and <^DE(i7) exists, the former will be dynamically selected in the dark 
energy dominated epoch. 17 We can therefore easily determine the trajectories from the 
existence conditions for A0MDE(a4) and A^>DE(a6). 18 In the strong coupling regime the 
existence conditions are very simple to lowest order; Ac/>MDE(a4) exists when Qm and Qa 
have similar signs, while Ac/>DE(a6) exists when Qa > 0. Thus each of the four trajectories 
occupies a region of parameter space of roughly equally large size. We shall comment on 
each of them later. First let us observe a common feature of all trajectories, namely that 
the vector field scales similarly as radiation, i.e., £Ia oc Q r = constant, during the radiation 
dominated epoch. This is due to the fact that the coupling term QaX on the right hand 
side of (2.47) is suppressed during RDE(i2). Hence the vector is effectively uncoupled and 
therefore scales as an ordinary Maxwell field, pa oc a -4 . 19 The fact that the vector is 
stable during the radiation dominated epoch suggests that inflation is a natural candidate 
for the initial conditions for the vector. Then, let us comment on each of the four types of 
trajectories. 

RDE(i2) —7- i^>MDE(i5) — > 0DE(i7) This trajectory is close to isotropy and charac- 
terize the standard coupled quintessence model [9]. See simulation a) of figure 2 for an 
example of this trajectory. It is realized if neither A^>MDE(a4) nor A^>DE(a6) exist, i.e., 
in the parameter region where 1 > ^QmQa and 2 > XQa- As explained above, the energy 
density of the vector field scales similarly as radiation during RDE(i2). During 0MDE(i5) 
and 0DE(i7) the vector typically decays rapidly, Q A oc a" 1 ^ 4 ^" 1 and Qa oc a ~^~ 2XQ ^ 
respectively. Since this type of trajectory exists in SCQ, the phenomenology at the per- 
turbative level is well understood [9]. 

RDE(i2) -> 0MDE(i5) -> A0DE(a6) This trajectory is close to isotropy in both the 
radiation and matter dominated epochs, but anisotropic in the dark energy dominated 
epoch. See simulation b) of figure 2 for an example of this trajectory. It is realized 
if A0MDE(a4) does not exist and A0DE(a6) exists, i.e., in the parameter region where 
1 £ 4QmQa and 2 < \Q A . During </>MDE(i5) the vector scales as £l A ~ a ~ 1+iQAQM and 
therefore decays rapidly for typical parameters consistent with the above condition. The 
trajectory will therefore typically have a transient period around the isotropic fix-point 
0DE(i7) (since some time is needed for Qa to grow) before settling at the anisotropic 
attractor A</>DE(a6). The trajectory is therefore usually more precisely described by 
RDE(i2) ->■ 0MDE(i5) -> ^DE(i7) ->■ A</>DE(a6). Notice that this is the case for sim- 
ulation b) of figure 2 where there is an isotropic period in the beginning of the dark energy 
dominated epoch. This period is short or entirely vanishing in the special case where 
^QmQa is just slightly smaller than unity. In general, however, the cosmology of this 

17 There may be a transient period close to </>DE(i7) before A</>DE(a6) is reached. An example of this is 
shown in simulation b) of figure 2. 

18 The isotropic fix-points are guaranteed to exists since we consider the case that both \Qm\ and A are 
smaller than unity. 

19 Here, we have assumed that the shear is negligible during the radiation dominated epoch. 
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a) Q_M=0.1, 1=0.1, Q_A=-10 




1/(10(1*10) 



b) Q_M=-0.01, A=0.1, Q_A=1000 




S/a6I 



c) Q_M=-0.01, A=0.1, Q_A=-1000 



3 




d) Q_M= 0.01, A=0.02, Q_A=1000 




E/a4E 

Figure 2. The phase flow and dynamics for four different sets of parameters (written in the figure), 
but with the same initial conditions close to the concordance radiation saddle RDE(i2). The four 
set of parameters give examples on each of the four different types of trajectories described in 
section 4.2. In each simulation (labeled a-d) the first figure in the corresponding column shows the 
phase flow (displaying the relevant fix-points that exists there), the second figure is a logarithmic 
plot of the time evolution of the dynamical variables, the third figure is similar with linear y-axis 
(it is then no point to plot the vector and shear since they are too small visually), and the fourth 
figure compares f2 m , Q r and fly to a ACDM simulation that starts with the same initial conditions 
(Si = \Qi/D,f CBM — 1| ). In the phase flow diagram the shear axis is normalized such that the fix- 
point that has the largest shear is positioned at x = 1 (for instance a4E is the shear at A(/)MDE(a4) 
which can be looked up in table 3). 
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trajectory is indistinguishable from SCQ up to today (fi^ ~ 0.7); the universe will first 
become significantly anisotropic in the future. 

RDE(i2) — > A0MDE(a4) — > 0DE(i7) This trajectory describes a universe which is 
close to isotropy during the radiation dominated epoch, anisotropic during the matter 
dominated epoch and close to isotropy again in the dark energy dominated epoch. In 
other words, anisotropization occurs during the transition to the matter dominated epoch, 
while isotropization occurs under the transition to the accelerated era. To the best of 
our knowledge, this is genuinely new behavior not predicted by any other known model. 
See simulation c) of figure 2 for an example of this. Such trajectories are realized when 
Ac/>MDE(a4) exists and A0DE(a6) does not exist, i.e., in the parameter region where 
1 < IQmQa and 2 > XQ A . 

RDE(i2) — > A0MDE(a4) — > A0DE(a6) This trajectory describes a universe which is 
anisotropic both in the matter and dark energy dominated epochs. Anisotropization occurs 
during the transition to the matter dominated epoch, while the transition to the dark 
energy dominated epoch represents a transition between two anisotropic scaling solutions. 
Also this is genuinely new behavior. Such trajectories are realized when both A(/)MDE(a4) 
and A0DE(a6) exists, i.e., in the parameter region where 1 < &QmQa and 2 < XQa- For 
an example of this trajectory, see simulation d) of figure 2. In that simulation the shear 
has approximately the same value during both A0MDE(a4) and Ac/)DE(a6). In general 
however, depending on the model parameters, the shear may be very different during the 
matter and dark energy dominated epochs. 

4.3 The present universe as a global attractor 

In the previous sections we have demonstrated that, given some conditions on the parame- 
ters and the initial condition, our model is, in many ways, dynamically very similar to the 
concordance model. We shall now see that also a more exotic, yet interesting, scenario can 
be realized. A global accelerated attractor with ~ 0.7 and O m ~ 0.3 is often discussed 
as an interesting possibility in the literature [79, 80]. In this scenario the present universe 
has already reached the global attractor and we do not live in the middle of a transition 
between two cosmic epochs. It is sometimes argued that such a scenario reduces the coin- 
cidence problem of the standard model, namely the question why we happen to live at a 
very special moment in the cosmic history [80]. 

The fix-point (i6)0 m is a global attractor with Q.^ ~ 0.7 and fl m ~ 0.3 when (Qm, A) = 
(—3.3, 2.2) and Qa < 0, see case b) of figure 1. In standard coupled quintessence it is 
impossible to realize a matter dominated epoch with these parameters since (/>MDE(i5) 
only exists for Q M < 3/2 [9]. In the frame work of that model one therefore needs to 
introduce a field dependent matter coupling, Qm(*) 5 which is small during <^>MDE(i5), but 
eventually grows large enough to realize the desired attractor [80]. In our model, however, 
there is no need to introduce time variation in Qm since the matter dominated epoch 
can be realized by A^MDE(a4) even with a large \Qm\- Note that A0MDE(a4) exists 
and has a small shear with Qm = —3.3 given the condition Qa <C Qm- In figure 3 we 
show how our model sucssessfully realizes the considered scenario through the sequence 
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RDE(i2) —7- A^MDE(a4) — > (i6)^ m . As far as we know there is no other model with a 
simple stationary coupling that can realize the present universe as a global attractor. 

Although our model can realize the present universe as a global attractor without any 
problems on the background level, it is expected that the large matter coupling will lead 
to problems for the perturbed universe. When the matter dominated epoch is realized by 
(/>MDE(i5) there is a strong bound \Qm\ < 0.13 which comes from a full CMB analysis [70]. 
As we shall see in section 4.5, it is expected that the bounds on \Qm\ is somewhat relaxed 
when the matter dominated era is realized by A^MDE(a4). Still we believe \Qm\ = 3.3 is 
way too large, but a full study of perturbations is required to verify this. 




Z/a4Z 



Figure 3. Simulation where the present universe is a a global attractor represented by (16)^ with 
~ 0.7 and Vl m ~ 0.3. The parameters used are Qm = —3.3, A = 2.2 and Qa = —100. 



4.4 Bounds from CMB quadrupole 

In this section we shall constrain the model parameters observationally using the CMB. 
Calculating the expectation values of the CMB anisotropies in the Bianchi I spacetime is 
in general a comprehensive task. Perturbation theory of the Bianchi I spacetime is to a 
large extent worked out (see [48, 81-83] and references therein), but a complete numerical 
implementation in a Boltzmann code is so far not achieved. When the shear is very 
small, however, one can treat the spacetime as a perturbation about the FLRW metric and 
consider only the leading corrections to the CMB. In this approach (used for instance in 
[48, 54, 56, 57, 84]) one consistently accounts for the anisotropic redshift of the photons 
emitted at the last scattering surface, but neglects any correction beyond this zero order 
effect. Note that the redshift in Bianchi models depends on the direction of propagation 
since the spacetime background is anisotropic. The total temperature anisotropy field 
AT/ (T) therefore consists of one part coming from inflation AT 1 / (T), which is similar 
as in the FLRW metric, and a "redshift part" AT"/ (T) due to the anisotropic metric, 
i.e., AT/ (T) = AT*/ (T) + AT a / (T). In appendix C we show that the leading term in 
AT a / (T) is a quadrupole term Q a which we shall refer to as the shear quadrupole. The 
quadrupole of AT*/ (T) we shall denote Q l and refer to as the inflationary quadrupole. All 
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shear multipoles higher than the quadrupole are strongly suppressed and can be neglected 
whenever the shear is sufficiently small to be compatible with observations. 

The total (observed) quadrupole Q depends on Q l , Q a and their mutual orientation. 
Depending on the orientation the total quadrupole may be larger or, more interestingly, 
smaller than the one coming from inflation. The 7-year WMAP data gives the value [29]: 

Q2 " 20 °(|^) 2 ~ (l0 ~ 5)2 ' (4.5) 
while the expectation value for the ACDM concordance model is: 

((QT^eg 2 . (4.6) 

Even taking into account the cosmic variance, y/2/5 ((Q*) 2 ), the observed value is supris- 
ingly low. It is well-known that for a shear quadrupole roughly of the same size as the 
inflationary quadrupole, the total quadrupole can be lowered to match the observed value 
[56, 57]. This requires a suitable orientation of Q 1 relative to Q a which, in fact, is not 
improbable, see [52]. 

Let us use this to constrain the parameters of our model. Although the possibility 
Q a ~ Q 1 is attractive, since it can explain the low observed quadrupole, it is clear that 
we must also allow parameters predicting Q a <C Q 1 . In the latter case the predicted total 
quadrupole is essentially equal to the one predicted by ACDM and the low observed value 
must be accepted as a statistical fluke. It is also clear that Q a cannot be several orders of 
magnitude larger than Q l since this would give, for any orientation, a large total quadrupole 
incompatible with observations. We shall use that the shear quadrupole must be around 
the same order of magnitude as the inflationary quadrupole or smaller: 

Q a < 10" 5 . (4.7) 

Then, let us relate the shear quadrupole to the parameters of our model. For con- 
venience we fix the coordinates such that the scale factors are unity today, i.e., a\\(to) = 
a±(to) = 1 and therefore a(to) = cr(to) = 0. The shear quadrupole is then given by (see 
appendix C for a derivation): 

S» = ^M. (4.8) 

where a& c = o"(idc) is the metric function defined in (2.17) evaluated at the decoupling 
time. Note that Q a is denoted Q2 in appendix C. To put bounds on the model parameters 
we shall derive an analytical expression for Q a . To derive a robust order of magnitude 
estimate we assume that the shear is given by fix-point A0MDE(a4) all the way from the 
decoupling time (z ~ 1100) up to today. 20 Since S is constant (as in any fix-point), it 
follows from definition (2.18) that a^c = Sa^c- 21 Inserting into (4.8), using the value for 



20 Since the transition to dark energy dominance must have happened quite recently, around z ~ 0.3, this 
is a quite good approximation although we know that the system (verified by simulations) will be closer to 
A0DE(a6) at z = 0. 

21 Note that, since |Qdc| — log 1100 is of order unity, |E| <C 1 implies |<Tdc| <S 1. 
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the shear in (3.4), we get 

_ 2|q dc | (AQ M Q A - 1) 

where the scale is related to the redshift by a^ c = — log(zdc + 1)- The constraint (4.7) then 
leads to: 

4gM ^ = 1 < ID" 5 , (4.10) 

where we have omitted prefactors. Note that this bound holds only if 4QmQa £ (lj°°) 
as required by the existence of A0MDE(a4). For typical parameters consistent with the 
existence of (a4) the constraint (4.10) can therefore be rewritten as: 



\Qm/Qa\ < nr 5 . (4.ii) 



For special parameters such that ^QmQa is close to unity one should use (4.10); otherwise 
(4.11) holds. 

So far we have neglected the shear during the accelerated era, £( a 6). Let denote the 
time of equivalence between matter and dark energy, i.e., £l$(t*) = O m (t*) ~ 0.5. Since 
t* corresponds to a low redshift, ~ 0.3, E(a4) will provide the main contribution to the 
shear quadrupole if E( a 4) ~ £( a 6)- If the matter era is realized by 0MDE(i5), however, like 
simulation b) of figure 2, or if X( a6 ) 3> £( a 4), then it is S( a 6) that will provide the main 
contribution to Q a . It is clear that we will get further bounds on the model parameters by 
considering the shear of A0DE(a6). Using the value for the shear in (3.7) as an estimate 
for the shear between t* and to we get the following estimate for the shear quadrupole: 



The constraint (4.7) then leads to: 



< ID" 4 , (4.13) 



Q 



.4 



This bound holds only if XQa £ (2, oo) as required by the existence of (a6). For typical 
parameters consistent with the existence of (a6) the constraint (4.10) can be rewritten as: 



|A/Qa| < 10- 4 . (4.14) 



For special parameters such that AQa/2 is close to unity one should use (4.13); other- 
wise (4.14) holds. Note that the bound on jA/Q^j is relaxed with an order of magnitude 
compared to \Qm/Qa\- 

4.5 Sound horizon 

In the previous section we showed that consistency with the observed CMB quadrupole 
requires |(2a/|/|Qa| ^5 10~ 5 if A0MDE(a4) is responsible for the matter era. The smallness 
°f |Qa/|/|Qa| have several consequences. Firstly it implies that the expansion is almost 
shear-free. Secondly it implies that the background expansion history, which goes as a tx 
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£(2/3)(i-Qm/6Qa) during A<^MDE(a4), is practically indistinguishable from the standard 
ACDM matter era where a oc t 2 / 3 . This is in contrast to 0MDE(i5) where a oc £ 2 /( 3+2< 2m) 
and the bound Qm < 0.13 [70] allows for a significantly different expansion history. 

The sound horizon at decoupling time is governed by the background expansion history. 
Assuming that the matter era is realized by a fix-point that includes the kinetic part of the 
scalar field, one can derive the following approximation for the sound horizon at decoupling 
[3, 9, 85]: 

r s ^r z^ MX , ( 4 - 15 ) 

where tq is the standard ACDM sound horizon, z^ c is the redshift at decoupling time and 

X is the dimensionless variable for the kinetic part of the scalar field defined in (2.41) 

-Q 2 

(which must be evaluated in the given fix-point). For 0MDE(i5) one gets r s = tqz a M . 
For Qm = 0.1, for instance, this gives a sound horizon 7% smaller than in ACDM. Since 
both photons and baryons are uncoupled in our model, we can calculate the sound horizon 
in exactly the same way. In particular, the exact background solution corresponding to any 
matter fix-point in our model, i.e., (3.1)-(3.3) with = /3 m , leads to the same expression 
for the sound horizon (4.15). For A0MDE(a4) (4.15) gives r s = roZ^ M ^^ A which means 
that the sound horizon is indistinguishable from the standard horizon. In the concordance 
model we know that the multipole location of the first acoustic peak is roughly inversely 
proportional to the sound horizon, / pca k oc r~ l . Conjecturing that the same relation also 
holds in our model, it follows that A^MDE(a4) cannot move the position of the first 
acoustic peak by means of a modified sound horizon. This is in contrast to ^>MDE(i5) 
which can move i max several percents. 

Since it is the fraction Qm/Qa that determines the background expansion history 
in A0MDE(a4), the magnitude of Qm itself is not dynamically relevant. It is therefore 
expected that the bounds on \Qm\ separately is somewhat relaxed compared to the bound 
\Qm\ < 0.13 in standard coupled quintessence. A full analysis of perturbations is required 
to confirm this. 



5 Summary and discussion 

In this work we have explored the cosmological consequences of a quintessence field with 
couplings to a Maxwell-type vector field as well as to dark matter. Both couplings have 
separately been studied extensively in the literature, but have until now not been considered 
simultaneously. In a dark energy context, motivated by higher-dimensional theories such 
as string theories, it is natural to consider the case where both types of couplings are 
present. We have showed that genuinely new behaviors arise due to the double coupling, 
and we have identified an interesting parameter region, corresponding to a strong vector 
coupling regime, yielding exciting and viable cosmologies close to the ACDM limit. Below 
we shall summarize our main results and discuss possible extensions and generalizations of 
the model. 

Doubly Coupled Quintessence (DCQ) can be viewed as a generalization of Standard 
Coupled Quintessence (SCQ) [9]. This is reflected in the dynamical system where all 
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the isotropic fix-points also exist in the simpler SCQ model. The dynamical system is 
characterized by 13 physically different fix-points, 7 isotropic and 6 anisotropic (see tables 1 
and 2 respectively). The properties of these fix-points are controlled by three dimensionless 
parameters A, Qa and Qm which determine the shape of the scalar potential and the 
strength of the couplings to the vector and matter fields, respectively. We identified a 
large parameter region, the strong vector coupling regime (see (4.1)), which yields viable 
cosmologies. In this region, the dynamics is characterized by only a handful of fix-points and 
can easily be understood qualitatively. In section 4.1 we established that viable cosmologies 
must have a radiation dominated epoch close to the isotropic concordance radiation saddle 
RDE(i2). This is required to have a sufficiently stable matter dominated epoch squeezed 
in between the radiation and dark energy dominated epochs. 

Interestingly, both the matter and the dark energy dominated epochs can be realized 
by anisotropic solutions, i.e., scaling solutions where a subdominant, but stable, vector field 
sources a small anisotropy in the expansion rate. A<^MDE(a4) is a genuinely new matter 
dominated scaling solution with a small stable vector field and a small quintessence field 
which can be responsible for the matter dominated epoch. A<^>DE(a6) is a quintessence 
dominated scaling solution with a small stable vector field. This fix-point was first discov- 
ered in the context of inflation [11, 14], but in the framework of our model it plays the role 
as the future attractor responsible for the late-time acceleration. There are also isotropic 
alternatives for the matter and dark energy dominated epochs, provided by <^>MDE(i5) and 
0DE(i7) respectively. These isotropic fix-points characterize the dynamics of SCQ. Thus 
we have anisotropic and isotropic alternatives both in the matter and the dark energy dom- 
inated epochs. Which trajectory that is dynamically preferred depends on the parameters 
in the theory, and we used the dynamical system approach to address this. In a certain 
parameter region where there are both anisotropic and isotropic alternatives for a given 
cosmic epoch, the anisotropic solution will be dynamically selected. It follows that the 
isotropic solutions are realized only in parameter regions where the anisotropic solution 
do not exist. When \Qm\ and A are smaller than unity there are four different types of 
trajectories which characterize the dynamics of the strong vector coupling regime: 

• isotropic both in the matter and the dark energy dominated epochs (equivalent to 
SCQ): 0MDE(i5) -> 0DE(i7), 

• isotropic matter dominated epoch and anisotropic dark energy dominated epoch: 
0MDE(i5) -»• A0DE(a6), 

• anisotropic matter dominated epoch and isotropic dark energy dominated epoch 
(isotropization occurs under transition to the dark energy dominated epoch): 
A0MDE(a4) -> 0DE(i7), 

• anisotropic in both the matter and dark energy dominated epochs: A0MDE(a4) — > 
A0DE(a6). 

All trajectories occupy roughly equally large regions of parameter space. Simulations for 
all types of trajectories are shown in figure 2, while further details are discussed in section 
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4.2. The rich dynamical behavior clearly distinguishes the model from other proposed 
late-time anisotropic models. In models where the anisotropic expansion is sourced by an 
isotropy- violating dark energy field [45-52], the anisotropization typically occurs during the 
transition to the dark energy dominated epoch (z < 0.3). This behavior is a reminiscent 
of the second special case </>MDE(i5) — > A</>DE(a6) of our model. Another possibility often 
considered in the literature is isotropy violation by a cosmic magnetic field (usually assumed 
to have a primordial origin) [53-60]. The energy density of such a magnetic field scales 
similarly as radiation (oc a -4 ) which means that the dimensionless density parameter will 
be stable during the radiation dominated epoch while it will decay (oc a -1 ) during the 
matter dominated epoch. Thus the anisotropy will die out during the matter dominated 
epoch and will be negligible at low redshifts. Although there are some similarities, this is 
clearly different from the third special case A</>MDE(a4) — > </>MDE(i5) of our model where 
the anisotropy during the matter dominated epoch arise during the transition from the 
radiation dominated epoch and then is stable all the way to the transition to the dark 
energy dominated epoch where it dies out. In the fourth type of trajectory the transition 
to the dark energy dominated epoch represents a transition between two anisotropic scaling 
solutions, also a new type of behavior. 

Let us also briefly comment on a fifth type of trajectory that occurs if we allow \Qm\ 
and A to be larger than unity. For (Qm, A) = (—3.3, 2.2) the isotropic fix-point (i6)^ m 
represents a global attractor with ~ 0.7 and Q m ~ 0.3. Thus we have the exotic 
possibility that the present state of the universe is the global attractor. This scenario is 
not possible to realize in the framework of SCQ since it is impossible to realize a matter 
dominated epoch. In our doubly coupled model, however, the matter dominated epoch 
will be realized by A</>MDE(a4). It is therefore possible to realize a cosmologically viable 
sequence with the present universe as the global attractor. Nevertheless, problems are 
expected at the perturbative level, see the discussion in section 4.3. 

The CMB quadrupole is very sensitive to shear and we used this to constrain the 
parameters of the model. More concretely we put upper bounds on the ratio of the strength 
of the matter coupling to the strength of the vector coupling, and on the ratio of the 
parameter in the quintessence potential and the strength of the vector coupling. As derived 
in 4.4 the upper bound is \Qm/Qa\ ^5 10 -5 if the matter dominated epoch is realized 
by A0MDE(a4), and [A/QaI ^ 10~ 4 if the dark energy dominated epoch is realized by 
Ac/>DE(a6). The surprisingly low quadrupole observed in WMAP is in fact not improbable 
if | Qm I Qa I or | A/ Qa I lays around the upper bounds such that the shear quadrupole is of 
the same order of magnitude as the inflationary quadrupole. Given a suitable orientation 
(which is not improbable) , the observed quadrupole will then be lower than the inflationary 
quadrupole in harmony with observations [52, 56, 57]. 

Our model is dynamically equivalent to ACDM in the limit Qm/Qa — > 0, X/Qa — > 0. 
In this limit the anisotropic solutions A0MDE(a4) and A<^DE(a6) become equivalent to 
the matter and dark energy dominated epochs of ACDM, respectively. Therefore, although 
DCQ is very different from ACDM qualitatively, the strong bounds on the model param- 
eters imply that the quantitative differences on the background level are extremely small. 
For example the mean scale factor evolves as a oc 

t (2/3)(i-Q M /QA) during A^MDE(a4) which 
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is practically indistinguishable from ACDM where a oc t 2 / 3 . As verified in the simulations 
in figure 2, there are only minor differences from ACDM on the background level for an 
epoch realized by an anisotropic solution (when the constraints on \Qm/Qa\ an d |A/Qa| 
are satisfied). A full study of perturbations is required for a more complete understand- 
ing of the phenomenological consequences of the anisotropic epochs. For 0MDE(i5) the 
phenomenology at the perturbative level has been explored in the literature and the con- 
sequences for structure formation and the CMB are well understood [9]. The upper bound 
on the matter coupling is \Qm\ < 0.13 given that 0MDE(i5) is responsible for the matter 
dominated epoch [70]. Interestingly, we expect that the bound on \Qm\ is somewhat relaxed 
compared to this if the matter dominated epoch is instead realized by A0MDE(a4) since 
in that case the sound horizon at decoupling time is indistinguishable from ACDM (see 
section 4.5). 

There are several natural ways to generalize the model. For instance it would be 
interesting to consider time variations in the coupling constants and/or other types of 
scalar potentials that exhibit tracking behavior [70, 86]. It would also be interesting to see 
if the model can be (minimally) modified to realize a non-standard radiation dominated 
epoch. Since the constraint on a primordial shear is much weaker than at late times [87], 
it is of particular interest to check if such a minor modification could allow for a radiation 
dominated era realized by one of the anisotropic scaling solutions (a2)^ 4( . r and (a3)^ mr . To 
better understand the observational signatures of DCQ a full analysis of the perturbations 
is necessary. This work will be carried out in the near future. 

Finally, it should again be emphasized that the Maxwell-type vector field considered in 
this paper represents a possible dark component of the universe and must not be confused 
with the photons of the standard model of particle physics. Note that if our considered field 
was identified with photons, the coupling to quintessence would correspond to an extremely 
rapid time variation in the fine structure constant, af s oc f((f>)~ 2 , clearly incompatible with 
observational constraints (the upper bound on time variation of the fine structure constant 
between decoupling time and today is at the one-percent level [88, 89]). This rules out the 
possibility that our considered vector field could be identified with a possible homogenous 
component of large-scale cosmic magnetic fields. 

To conclude, we have showed that the model is dynamically very rich, and can realize 
various scenarios that violate isotropy while being close enough to the ACDM limit to 
be taken seriously. An important lesson learned is that although the universe looks like 
ACDM based on present observations, it certainly does not mean that the universe is 
ACDM. Cosmological observations have forced us to introduce two new dark fields beyond 
the standard model of particle physics. The present energy density of both dark energy and 
dark matter exceed the energy density of baryons by approximately one order of magnitude. 
We should therefore not be surprised if future observations lead to discoveries of other less 
dominant types of fields with too little energy or too weak couplings to be seen in present 
data. 
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A Eigenvalues 



The stability analysis in section 3 is based on the following eigenvalues of the matrix of 
linear perturbations around the fix-points: 

Fixpoint (il±) : (2, 0, 2 T 2v / 6Q A , 3 ± V6Qm, 6 ± V6X). (A.l) 
Fixpoint (i2): (4, -1, -1, 1, 0). (A.2) 



Fixpoint (i3): »|. (A.3) 

Qm 2 2Qm 2 2Qm Qm 



P«pomt(i4): (^-1,8 X , , ,-- + ^ J. (A.4) 

Fixpoint (i5): 

( ~ \ + Qm, ~\ + Qm, "I + 2Qm, -1 + 2Q m (2Qa + Qm), 3 + 2Q M (Q M - \)j . (A.5) 
Fixpoint (i6): 

( 3 | 3Q A/ l | 3Q M _ 1+ 3(2Q a + Q m ) 



2 2(A — Qm) A — Qm' A — Qm 

3(A - 2Q M ) 



4(A - Q M ) 

Fixpoint (i7): 



/ 8(3 + \{Q M - A))(3 + 2Qm(Q m - A)) 
V 3(A-2Q M )2 



1 - i X 2 - 2Q A X, 3 + Q M X, 6 + > 



Fixpoint (a2): - -i ± ±y/A±BC) , 

where: 

A = (32 + 48Qj - 8Q4 (2 + SQl) A - (7 + 12Q|) A 2 ) , 



j 



B = -4J{2 + 3Ql), 



(A.6) 



\2 \2 

-3 - XQ M + A 2 , -3 + — , -3 + y , -4 + A 2 , -4 + 2AQ A + A 2 ) . (A.7) 



(A.9) 



C = Jl2Q\\ 2 - 8Q A X (-4 + A 2 ) + 2 (-4 + A 2 ) 2 + 12Q3 A (4 + A 2 ) + 3Q 2 (16 + A 4 ). 
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Fixpoint (a3): U - -\± -±-VA±b) , (A. 10) 

V Qm 2 4Q m J 



where: 

A = 4 + 6Ql - 8Q A (2 + 3Q 2 A ) Q m - 4 (l + 3Q|) Q^, 



B = 2^ (2 + 3Q2) ( 2 + 3Qj + 8 g A (2 + 3Qjf) Q M + 4Q| f (_ 2 + 9Qj + i 2 Qj) + C) , 

C = 16Q A (-2 + 3Ql) Q|f + 4 (2 + 3Ql) Q^. 
Fixpoint (a4): 

^3 - 3A (6 + 9Ql + 7Q a <9a/) , 8 - 6A (6 + 9Qi + 7Q A Q M ) , 

^(-1 + Qm(Qa + 3Q m )A±ABC), (A.ll) 

3A (4 + 6Ql + 6Q a Qm + 4Qm - QaA - 3Q M A) , 
where: 

A = 1/ (4 + (2Q A + Qm)(3Qa + Qa/)) , 



B = -2 J (-2 - 3Qi - 2Q a Qm + Q 2 M ) , 



C = ^J{-8 + 32Q\Q M + 13Q M + 2QaQm (9 + + (-11 + 32Q&)) . 

Fixpoint (a5): extremely lengthy algebraic expressions in the parameters (Qm,Qa, A). 

Fixpoint (a6): 

^ - 3,4(8 + (6Q A - X)(2Q A + A)), -8A (4 + 6Ql + Q A A - A 2 ) , 

- - + 3A\(2Q A + A)8 + (2Q A + \)(6Qa + A) ± AS, (A.12) 

- 3,4(8 + (6Q A + 4Q M - 3A)(2Q A + A)) j , 
where: 

yl = l/(8 + (2Q A + A)(6Q A + A)), 

B = -^(-8 - 12Q| _ 4 q a a + A 2 ) (-72 + (2Q A + A)(17A + 2Q A (-19 + 4A(2Q A + A)))). 
B Non-viability of fix-point (a5) 

In section 3.3 we stated that (a5) cannot be responsible for neither the matter dominated 
nor the accelerated epoch. Here we show why. 

First we consider the case A ^> Qm in which case uj c s — 0. Let us check if this 
matter dominated scenario is cosmologically viable. In the limit A S> Qm the shear is 
E ~ — j + 1 3^-. Since the eigenvalues of the Jacobian matrix are extremely complicated, it 
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is convenient to neglect higher order terms in the small quantities a = b = — j + | 
Note that o <C 1 is required to have uj c g ~ 0, while b <C 1 is required to have a small shear 
|E| <C 1. To first order in a and 6 the conditions for existence are: 

a + b>0, (B.l) 
-1 + ^-A 2 - -^A 2 (2a + 6) >0, (B.2) 

2 - A 2 (a - -b) > 0. (B.3) 

When A is of order unity, this simplifies to A 2 > 3 and a + b > 0. For any A 2 > 3 it is 
possible to satisfy (B.1-B.3) simultaneously for small a and b. Thus we have showed that 
it is possible for (a5) to represent a matter dominated solution with a small shear. To 
say whether this possibility is cosmologically viable we need to know the stability. In our 
considered approximation the real part of the eigenvalues are: 

■l + 3a, -^ + ya + 46, -4(a + 6), ~(l-a), -|(1 - a)J (B.4) 

Here we have assumed A 2 > %r to simplify the eigenvalues (several terms then become 
imaginary). The existence conditions allows for a slightly smaller value, A 2 > 3, but from 
O m = 1 — A + 0(a) + 0(6) one see that A 2 S> 3 is needed to have matter as the dominant 
component. 22 Notice that the third eigenvalue is negative when (B.l) is satisfied. The other 
eigenvalues are also negative in our considered approximation. It follows that the matter 
dominated solution with a small shear is a stable attractor. The solution can therefore 
not escape from the matter era into the accelerated era. Thus (a5) does not provide a 
cosmologically viable alternative to the matter dominated era. 

Next we check if the fix-point can drive acceleration. Acceleration is realized if 

Waff = -7^-6 [-1,-1/3). (B.5) 
Qm — A 

Since we consistently assume A > in this paper, we have acceleration if and only if 
Qm < — 5 A. 23 It is possible to satisfy the conditions for existence together with the 
condition for acceleration only if Qm < —y/2/2. Thus a small matter coupling, \Qm\ <C 1 
is incompatible with (a5) driving acceleration. For the case Qm < -V2/2 there is a 
restricted region of parameter space that is compatible with the conditions for existence 
and acceleration. As we shall now see, however, this region does not provide a viable 
cosmology. 

First we check the interesting special case of an "almost de Sitter" solution, i.e., oj c q ~ 
— 1 and |S| <C 1. It is convenient to use the three quantities Qm, a = —1 — f^j and 
b = We assume that the latter two are small (a <S 1, b 1) which is equivalent to 
E <C 1 and u) e g ~ — 1. From the constraint O m > it follows that —2 — Q M {^a + b) > to 



22 For A 2 ~ 3 the dominating component is the scalar field with equal kinetic and potential part such that 
the total equation of state is zero. 



23 Note that for Qm > A we have u> c g < —1 and the fix-point do not exist. 
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first order in a and b. Thus the solution requires a large matter coupling, Q?j 3> 1. This 
rules out the possibility for (a5) to represent an "almost" deSitter solution. 

Next we consider the possibility of an accelerated scaling solution with ^ and 
7^ 0- I R this case it is convenient to use the three quantities Qm, ^cS an d X. Recall 
that we consistently consider the case A > 0. In terms of our considered quantities this 
translates into Qm < 0. Let us assume acceleration and a small shear, i.e., u! c g G (— |, — 1] 
and X <C 1. Given these assumptions the constraint > is automatically satisfied while 

2 

the conditions > and fl m > are equivalent to X > and Q\ t > 3 1+ ^f ■ The latter 
condition implies a minimum value on \Qm\- Tuning the model to the concordance model 
by setting uj c g = —0.7 implies > 4.9. In order to have a small shear we must require: 



A - GQa - 4Qm 
4(Qm - A) 



< 1. (B.6) 



Using io c ff = —0.7 it follows from (B.6) that Qm/Qa — —42/31 which is incompatible 
with a matter dominated epoch via A^MDE(a4). Also note that the condition > 4.9 
derived above is incompatible with a matter dominated epoch via ^MDE(i5) (the existence 
condition for 0MDE(i5) is Qfij < 3/2 ). Thus it is impossible to realize a matter dominated 
epoch in the parameter region where (a5) represents an accelerated epoch. This rules out 
(a5) as a candidate for the matter dominated era. 

C CMB quadrupole in the Bianchi type I spacetime 

Here we shall derive the leading effects of the shear on the CMB. Other derivations can be 
found in [48, 63]. Let us consider the most general Bianchi type I spacetime: 

ds 2 = -dt 2 + a{t) 2 dx + b(t) 2 dy 2 + c(t) 2 dz 2 . (C.l) 

Our first goal is to derive a formula for the redshift of photons coming from the last 
scattering surface. Since the expansion is anisotropic the redshift will be a function of the 
direction of the incoming photon: 

z(p) = ^, (C.2) 

*dc 

where \d c and Ao is the wavelength of photons at decoupling (last scattering surface) and 
today, respectively. We will derive this formula using an analogy to photons propagating 
in the flat FLRW metric. First consider the geodesic equation: 

^ = -I>V, (C3) 
where the only non-vanishing Christoffel coefficients are 

r° = ajdj, = r* = . (C4) 

Note that if the tangent vector of a photon is parallel to a coordinate line initially, the 
photon will continue to propagate along that coordinate line. For example if the photon's 
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four- velocity initially is parallel to the x-axis, u 2 = u 3 = 0, then = ^jx = according 
to the geodesic equation. We shall now use this observation to find an expression for the 
redshift as a function of the direction of the received photon. Let the photon's direction of 
propagation be defined by the unit (three) vector 

P % = iPxiPy,Pz) = (sin 6 cos <fr, sin 6 sin (f>, cos 8). (C.5) 

We now introduce new coordinates x^ = (t,x,y,z) which are (spatially) rotated relative 
to such that the x - axis is parallel to the photon's direction of propagation: 



x 



PxX+PyV+PzZ. (C.6) 



We let a, b and c be the scale factors in the new coordinate chart. The coordinate trans- 
formation defined by (C.6) determines the metric component g 11 : 

% 1 

According to the discussion above the photon will propagate along the x coordinate line. 
The photon's propagation is therefore subject to the null-condition: 

ds 2 = -dt 2 + d 2 dx 2 . (C.8) 

This is mathematically equivalent to a photon propagating in the flat FLRW metric. Thus 
the redshift is 

l + z(p) = ^, (C.9) 
a 

where a is the scale factor at the time of emission and do today. For convenience we use 
flo = = c o = 1 which implies do = 1. From (C.7) and (C.9) we then get an expression for 
the redshift as a function of the direction of propagation for the incoming photon emitted 
at the decoupling time: 

i+ ^/?S' (c - io) 

where the subscript indicates that the scale factors are evaluated at the decoupling time. 
This is our desired formula, the redshift parametrized by the direction of the photon. From 
here we follow more or less the derivation in [48]. It is convenient to rewrite the redshift: 

l + z(p) = ±y/l+p*e+pl7i, (C.ll) 

where 

.5) "=0 - 1 < C12 > 

Note that |e| and \rj\ essentially are the "eccentricities" squared. For a small shear the 
eccentricities will be small and we can expand to first order in e and rj. First we calculate 
the mean redshift over all directions to first order in e and rj: 

<*(p)> = ^ / dn p z(p) = -1 + i + 0{e 2 , r/ 2 ), (C.13) 
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where a = (06c) 1 ' 3 is the geometric mean of the scale factors. In the main text we defined 
an isotropic redshift z, see (2.22). Notice that to first order in the perturbations (z(p)) = z. 
Thus our defined isotropic redshift has an interpretation as the averaged redshift when the 
shear is small. 

The homogenous background temperature field will therefore be anisotropic: 

Tdc 



T a (p) 



l + *(p) (C . 14) 
aT dc (l-^ple-^p 2 zri ) + 0(e 2 , V 2 ), 



where T^ c is the decoupling temperature. In the following we shall omit the notation 
0(e 2 ,i] 2 ), but it is clear that we will neglect all second order terms. The total temperature 
field is T = T a + AT* where AT* is the (inhomogeneous) part coming from inflation. The 
average (total) temperature is: 

{T) = hj d °p Ta (p) = aTdc ( x - V e + ^) ■ (°- 15 ) 

The shear temperature anisotropy field is therefore: 

ms ^ = 1 -Tffi = lp + lfo-^ +v) , (c , 6) 

or in terms of the angular variables: 

6(0,6) = - e sin 2 8 sin 2 6 + -rj cos 2 9- -(e + r>). (C.17) 
2 2 6 

We shall now expand the field in the spherical harmonics: 

1 

5(9,6) =J2 ai m Y lm (8,6), (C.18) 

l m=—l 

where, due to the orthogonality of the spherical harmonics, the coefficients are given by: 



a lm = I dn p 6(p)Yi* m - (C19) 

The multipoles are defined 



Qt 



1 i(i + i) ' 



\ 2vr (21 + _ 

\ m=—l 



To first order in e and 77 there are only non- vanishing coefficients for 1 = 2: 



1 / 7T TT 

«20 = oA/t( 2? ? ~ e )' 021 = 02-1 = 0, a 22 = a 2 -2 = -W— e, (C.21) 
6 V 5 V oO 



which gives 



gi = ^- 7 =V / e 2 + ^ 2 -^- (C22) 
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Thus for our axisymmetric spacetime (2.17) the quadrupole is 



21 



Qi 



12 



(C.23) 



5^3 



0"dc 



To second order in our small quantities we also find a hexadecapole: 




(C.24) 



Note that Q% ~ {Qlf ■ More generally, for even multipoles we have Q% n ~ (Q2)" while 



anisotropy field is a quadrupole coming from the shear that will mix with the inflationary 
quadrupole Q\. Since observationally we know that Q2 < 10 -5 it is clear that we can 
neglect the correction to the hexadecapole and any higher multipoles. 
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